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In  the  first  phase  of  this  project,  wc  constructed  experimental  prototypes  of  the  contiguous 
domain  oscillator  device,  'fhese  prototype  devices  oscillated  in  the  single-domain  mode  (similar 
to  a  lateral  Gunn  diode),  producing  signals  in  the  range  from  6  to  28  GHz.  Contiguous  domain 
operation  could  not  be  achieved.  We  believe  this  failure  was  due  to  the  high  sheet  resistance  of  , 
our  resistive  gate.  (With  the  gate  resistance  high,  image  charge  motion  in  the  gate  disrupts  the  ' 
ideal  linear  gate  potential). 

The  objective  of  the  second  pha.se  of  research  was  to  produce  new  devices  having  lower  gate 
resistance  and  to  demonstrate  contiguous-domain  operation,  Tliese  new  devices  were  to  be  built 
by  nr  Gallium  Arsenide  Technology  Center  and  were  to  be  delivered  to  us  by  raid- 1991.  Due 
to  a  scries  of  problems  at  ITT,  the  devices  were  never  delivered.  Instead,  we  have  shifted  our 
emphasis  to  computer  simulation,  and  have  developed  a  new  computer  model  of  the  device 
which  shows  the  transition  from  multiple  to  single-domain  mode  as  the  gate  is  moved  further 
away  from  the  channel.  This  new  simulation  program  allows  us  to  study  the  transition  legime 
and  to  investigate  the  range  of  device  parameters  over  which  contiguous-domain  operation  can 
be  expected. 
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Background 

llic  contiguous  domain  oscillator  (CDO)  is  a  resistive-gate  GaAs  FET  which 
functions  as  a  \oltagc-controllcd  millimctcr-wave  oscillator,  llie  structure  is 
compatible  with  planar  processing  and  can  be  incorporated  with  standard  GaAs 
FT.Ts  in  the  fonn  cf  monolithic  millimeter-wave  integrated  circuits  (MiMICs).  17ic 
device  was  first  proposed  by  Cooper  and  Tlmmbcr  in  1983,  and  has  been  under 
investigation  at  Purdue  since  1985. 

Conventional  microwave  oscillators  (such  as  the  IMPATT,  BARITT,  and  Gunn 
diodes)  generate  microwave  power  by  using  a  mechanism  which  depends  strongly 
on  electric  field  (e.g.  avalanche  multiplication  or  negative  differential  mobility)  to 
create  charge  packets  or  domains.  ITie  conventional  devices  all  employ  a  two- 
terminal  geometry,  as  shown  in  Fig.  1,  and  hence  are  diodes.  Since  tlie  internal 
electrostatic  geometry  is  one  dimensional,  once  a  charge  packet  is  introduced,  it 
alters  the  electric  field  throughout  the  device  (Fig.  la)  in  such  a  way  as  to  turn  off 
the  generation  mechanism.  The  single  charge  packet  tlien  drifts  across  the  device 
and  is  extracted  by  the  anode.  Once  the  charge  packet  has  drifted  out  of  the  device, 
the  internal  field  returns  to  it’s  original  value  and  the  generation  process  begins 
again.  The  entire  process  is  therefore  a  series  of  repeated  transients,  with  the 
oscillation  frequency  determined  by  the  inverse  of  the  generation  time  plus  the 
transit  time.  Clearly,  once  the  device  has  been  built,  the  transit  time  (and  hence  tlie 
oscillation  frequency)  is  fixed.  In  addition,  for  millimeter-wave  operation  the  drift 
distance  must  be  made  very  short,  typically  less  than  1  pm. 

In  contrast,  the  contiguous  domain  oscillator  employs  a  two-dimensional 
electrostatic  geometry,  as  illustrated  in  Fig.  lb.  In  this  stracture,  the  electric  field 
within  the  device  is  created  by  the  potential  drop  along  a  resistive  "gate"  electrode, 
and  electrons  are  constrained  to  drift  parallel  to  the  resistive  gate  by  a  potential 
barr’cr,  shown  schematically  in  Fig.  lb.  A  charge  packet  or  domain  in  this 
structure  alters  the  electric  field  only  locally  -  outside  the  packet  the  field  is  still 
determined  by  the  potential  gradient  in  the  resistive  gate.  If  electrons  are  supplied 
by  a  source  contact  and  extracted  by  a  drain  contact,  the  structure  in  Fig.  lb 
resembles  a  field-effect  transistor  with  a  resistive  gate.  Computer  simulations  [1,2] 
indicate  that  when  the  electric  field  is  in  the  regime  of  negative  differential 
mobility,  charge  domains  will  form  spontaneously  in  the  channel  near  the  source 
and  drift  along  the  channel  into  the  drain.  In  contrast  to  microwave  diodes,  the 
frequency  is  not  determined  by  the  physical  length  of  the  device,  but  rather  by  the 
spatial  period  of  the  domains  in  the  channel.  This  spatial  period  is  found  to  depend 
upon  the  average  electron  density  in  the  channel,  which  can  be  controlled  by  the 
source  to  gate  voltage. 


E{X) 


(a) 

Gunn 

Diode 


Figure  1.  Schematic  illustration  of  the  operation  of  (a)  the  conventional  Gunn 
d.ode  (top),  and  (b)  the  Contiguous  Domain  Oscillator  (bottom).  The  Gunn  diode 
l.as  a  one-dimensional  electrostatic  geometry,  and  the  presence  of  a  chrage  packet 
perturbs  the  electric  field  throoughout  the  device.  The  CDO  has  a 
two-dimensional  electrostatic  geometry,  and  the  charge  packet  does  not  perturb 
the  electric  field  in  the  remainder  of  the  device. 
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Experimental  devices  were  fabricated  both  at  Purdue  and  by  the  ITT  Gallium 
Arsenide  Technology  Center  (ITT-GTC)  in  Roanoke,  VA.  A  schematic  cross 
‘Section  of  the  experimental  contiguous  domain  oscillator  device  is  shown  in  Fig.  2. 
h''  early  devices  fabricated  at  Purdue  failed  to  operate,  but  devices  delivered  to 
Purdue  by  ITT-GTC  in  February  1988  eventually  produced  microwave  signals  in 
the  frequency  range  between  6  and  28  GHz.  Although  these  frequencies  reach 
almost  to  the  millimeter-wave  range,  they  are  too  low  to  be  due  to  contiguous 
domain  oscillation.  We  were  able  to  show  that  the  device  was  operating  in  the 
single  domain  mode:  in  effect,  only  one  domain  was  forming  in  the  channel  at  a 
time,  and  this  domain  propagates  from  source  to  drain  before  a  new  domain  can 
form.  The  frequency  is  simply  the  inverse  of  the  transit  time  from  source  to  drain, 
and  scales  inversely  with  channel  length,  as  shown  in  Fig.  3.  Some  frequency 
modulation  (up  to  20%)  could  be  obtained  by  adjusting  the  gate-to-source  voltage, 
but  not  as  much  as  expected  from  the  contiguous-domain  mode. 

We  suspect.. .  that  the  problem  in  the  initial  devices  was  due  to  the  relatively  high 
resistivity  of  the  resistive  gate.  This  could  be  demonstrated  by  a  simple 
calculation:  We  first  assume  an  ideal  gate,  in  which  the  potential  increases  linearly 
from  source  to  drain.  The  operation  of  this  ideal  device  is  simulated,  and 
contiguous  domains  are  observed  in  the  channel,  as  expected.  We  then  assume  that 
the  charge  domains  in  the  channel  induce  equal  and  opposite  image  charges  on  the 
resistive  gate,  and  that  tliese  image  charges  move  along  with  the  domains.  The 
charge  motion  in  the  gate  constitutes  a  current  flow,  and  Jiis  current  multiplied  by 
the  gate  resistivity  results  m  local  voltage  drops.  When  these  local  potential  drops 
are  added  to  the  ideal  linear  potential  drop,  the  desired  linear  voltage  drop  is 
severely  perturbed  -  the  lateral  electric  field  in  the  gate,  instead  of  being  constant, 
oscillates  sinusoidally  with  position,  and  the  minimum  field  actually  becomes 
negative  at  the  low  point  of  each  cycle!  Such  a  severe  perturbation  must  destroy 
the  screening  effect  of  the  gate: 

The  obvious  solution  was  to  decrease  the  resistivity  of  the  resistive  gate.  The 
original  ITT  devices  have  a  WSiN  gate  with  a  sheet  resistivity  of  30  kfi  per  square. 
Our  calculations  indicated  that  this  needs  to  be  reduced  to  about  3  per  square  to 
eliminate  the  local  field  perturbations.  It  then  became  the  primary  objective  of  the 
research  to  build  new  devices  meeting  these  specifications 

ITT  initially  agreed  to  supply  six  three-inch  wafers  of  CDO  devices  by  mid-1991. 
All  were  to  have  gate  resistivity  between  1.5  and  3  kO  per  square.  However,  due 
to  a  series  of  internal  problems  at  ITT,  beginning  in  the  summer  of  1991,  the 
wafers  were  delayed.  These  problems  were  both  technical  and  (apparently) 
political.  In  spite  of  numerous  attempts  to  stimulate  delivery  of  the  samples  via 
phone  calls,  FAX  messages,  etc.,  we  were  eventually  forced  to  conclude  that  the 
devices  were  not  going  to  be  forthcoming. 
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Resistive  Gate:  WSiN,  110  nm,  30  kn/  □ 
N  Channel;  1.2x10^’'cm'^.  150  nm 
S  and  D:  SxIO^^cm'^  150  nm 

Distance  X :  {  q 

Figure  2.  Cross  section  of  the  experimental  devi( 
structure  is  essentially  a  resistive-gate  MESFET, 

GaAs  rex  processing. 


Figure  3.  Measaured  oscillation  frequency  plotted  vs  inverse  length. 
The  linear  relationship  indicates  that  the  frequency  is  the  inverse  of  the 
transit  time,  and  proves  that  the  devices  are  operating  in  the 


single-domain  mode. 
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As  a  result  of  ITT’s  de-facto  withdrawal  from  the  project,  we  redirected  the 
program  at  Purdue  to  concentrate  on  obtaining  a  more  thorough  understanding  of 
the  operation  of  the  device  under  circumstances  where  both  the  source/drain 
electric  field  and  the  resistive-gate  electric  field  exert  competing  influences  on 
electron  motion  in  the  channel. 


Funding  History 

The  original  CDO  devices  were  built  and  tested  under  grant  no.  AFOSR-85-0193, 
which  ran  from  1  May  1985  through  31  April  1988.  A  continuation  grant 
(AFOSR-85-0193E)  and  it's  no-cost  extension  supported  one  graduate  student  from 
1  May  1988  through  31  January  1990.  From  1  Februaiy  1990  through  30  June 
1991  we  conducted  this  research  using  local  money,  primarily  a  grant  from  the 
Indiana  Corporation  for  Science  and  Technology.  With  these  local  funds  depleted, 
we  applied  for  and  received  the  present  one-year  grant,  AFOSR-9 1-0224.  This 
grant  began  on  1  July  1991  and  terminated  on  30  June,  1992.  A  no-cost  extension 
to  31  December  1992  was  requested  and  approved. 


Device  Simulation  Results 

In  the  original  computer  program  used  to  simulate  the  contiguous  domain  effect, 
the  local  potential  at  each  point  in  the  channel  is  calculated  using  one-dimensional 
electrostatics.  The  potential  then  depends  on  only  two  variables:  the  local  gate 
voltage  and  the  electron  density  in  the  channel  at  that  point.  (A  non-local 
correction  term  is  also  included  to  account  for  fields  from  adjacent  regions  having 
drastically  different  electron  densities  —  this  eliminates  the  so-called  "gradual 
channel  approximation".) 

The  original  approach  is  reasonable  under  the  assumption  that  the  channel  potential 
is  dominated  by  the  resistive  gate,  and  indeed  has  been  verified  experimentally  for 
resistive  gate  devices  in  silicon  (work  performed  by  the  author  at  Bell  Labs  prior  to 
1983).  However,  this  one-dimensional  approach  is  not  valid  in  situations  where 
the  local  potential  is  influenced  by  the  potentials  on  the  source  and  drain,  i.e.  where 
the  gate  does  not  dominate.  It  is  precisely  these  conditions  which  give  rise  to 
single-domain  operation.  Thus,  in  order  to  investigate  single-domain  operation  and 
the  transition  between  single-  and  multiple-domain  behavior,  it  was  necessary  to 
write  a  completely  new  simulation  program  using  a  different  approach. 
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The  difference  between  the  bid  and  new  programs  is  illustrated  by  Figs.  4  and  5. 
In  both  circuits,  the  Jj  elements  at  each  node  i  represent  the  current  due  to  electron 
motion  in  the  channel,  'fhis  current  is  calculated  as  follows.  First,  the  local 
channel  potential  (|)i  is  calculated  at  each  grid  node  (more  about  this  later). 
Knowing  the  local  potentials,  the  electric  field  between  each  node  is  determined. 
The  currents  Ji  are  then  calculated  based  on  the  nonlinear  velocity-field 
relationship  for  electrons  in  GaAs.  A  time  step  At  is  taken,  and  the  local  electron 
densities  are  allowed  to  change  as  a  result  of  the  currents  flowing  into  and  out  of 
each  node.  The  new  electron  densities  a'e  then  used  to  calculate  new  potentials, 
and  the  cycle  repeats.  If  the  discretization  in  space  and  time  are  sufficiently  small, 
the  difference  equations  are  good  approximations  to  the  differential  equations 
governing  device  operation,  and  valid  results  are  obtained. 

The  difference  between  the  old  program  (Fig.  4)  and  the  new  program  (Fig.  5)  is 
the  manner  in  which  the  local  potentials  <j)i  are  calculated.  Although  the  figures 
look  very  similar,  the  only  difference  being  the  inclusion  of  the  lateral 
semiconductor  capacitances  Cs  in  Fig.  5,  this  difference  is  extremely  important.  In 
fact,  tlie  approaches  used  to  calculate  the  potentials  for  the  two  cases  are  totally 
different.  In  the  original  program  (Fig.  4),  the  local  potentials  (j)i  could  be 
calculated  from  the  local  gate  potentials  Voi  and  the  local  electron  densities  ni 
through  a  simple  algebraic  equation.  In  the  circuit  of  Fig.  5,  the  local  channel 
potentials  ({)i  must  be  calculated  by  solving  the  set  of  simultaneous  equations 
describing  the  entire  device.  This  requires  the  use  of  a  simultaneous  linear 
equation  solver  at  each  time  step.  Since  these  equations  include  the  source  and 
drain  potentials  implicitly,  the  source  and  drain  may  dominate  over  the  gate  if  the 
element  values  in  the  model  are  of  appropriate  values. 

In  the  model  of  Fig.  5,  if  the  coupling  between  the  gate  and  channel  is  reduced,  for 
example  by  moving  the  gate  further  away  from  the  channel,  the  Cd  capacitors 
connecting  the  gate  and  channel  become  small  and  the  channel  potential  will  be 
dominated  by  the  source  and  drain  potentials.  The  device  then  functions  as  a 
lateral  Gunn  diode,  and  the  model  produces  single* domain  oscillation.  On  the 
other  hand,  if  the  gate  is  placed  in  close  proximity  to  Lhe  channel,  as  in  the  CDO, 
the  gate  dominates  and  the  model  of  Fig.  5  produces  the  same  results  as  the 
original  model  of  Fig.  4. 

Figure  6  shows  the  electron  density  as  a  function  of  position  for  several  times  for 
the  circuit  of  Fig.  5  when  the  gate  has  been  placed  at  a  large  distance  (600,000  A) 
from  the  charmel.  Here  the  device  is  functioning  essentially  as  a  lateral  Gunn 
diode,  with  the  electric  field  lines  from  the  charge  domains  terminating  on  the 
source  and  drain,  and  a  single  dipole  domain  propagates  along  the  channel  toward 
the  drain.  Multiple  domains  do  not  form. 


Figure  4.  Equivalent  circuit  illustrating  the  original  program  for  simulation  of  the 
contiguous-domain  effect.  The  full  simulation  consisted  of  256  nodes,  not  just  the 
five  shown  here.  The  channel  potential  0  is  calculated  at  each  node,  based  only 
on  the  gate  potential  and  electron  density  at  that  node.  Once  the  channel 
potentials  are  known,  the  currents  J  are  calculated  and  a  time  step  taken. 


Figure  5.  Equivalent  circuit  illustrating  the  new  simulation  program.  Here  the 
chf»nnel  potentials  are  obtained  by  simultaneous  solution  of  the  node  equations  of 
the  entire  circuit.  In  this  way,  the  effect  of  source  and  drain  potentials  are 
included  in  a  natural  way.  If  the  coupling  capacitance  Cd  between  gate  and  drain 
are  small  compared  to  the  semiconductor  capacitances  Cs  connecting  each  Ay 
region  of  the  device,  then  the  source  and  drain  field  will  dominate  and  single 
domain  oscillations  will  be  observed. 


‘jgl=0  ^  vgg=8^  vs=0 . 5^  vd=8 . 5^  iox=600000A 

Figure  6.  Electron  density  in  the  channel  as  a  function  of  position  predicted  by 
the  new  simulation  program  when  the  gate  is  moved  well  away  from  the  channel 
so  that  the  source/drain  field  dominates.  A  single  dipole  domain  is  established  as 
an  initial  condition  at  t=0.  This  domain  stabilizes  and  drifts  along  the  channel,  as 
shown  at  t=10,  40  and  80  ps.  This  type  of  operation  is  referred  to  as  the  SINGLE 
DOMAIN  MODE. 
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In  contrast,  Fig  ows  the  predictions  of  the  same  program  when  the  gate  is 
brought  into  ck  ^  jximity  to  the  channel  (1050  A).  Here  we  plot  electron 
density  for  one  particular  time,  150  ps  after  the  start  of  the  simulation.  The 
electron  distribution  in  the  channel  has  formed  multiple  domains  in  exactly  the 
manner  predicted  by  our  original  program.  The  oscillation  frequency  observed  in 
the  drain  current  for  this  case  is  63.2  GHz. 

The  new  simulation  program  has  been  compared  with  the  original  program  under 
conditions  where  the  gate  dominates  the  channel  potential,  and  the  predicted 
oscillation  frequencies  are  in  very  good  agreement.  TTiis  is  shown  in  Fig.  S,  where 
frequency  is  plotted  as  a  function  of  the  product  of  electron  density  ns  and 
depletion  depth  xd-  The  small  symbols  are  predictions  of  the  original  program  for 
a  large  numl^r  of  simulation  runs,  and  the  large  black  circles  are  predictions  of  the 
new  program  for  several  runs.  In  spite  of  differences  in  the  models  and  the 
algorithms  used  to  calculate  them,  the  predictions  of  the  two  approaches  are  in 
4uite  good  agreement.  This  agreement  gives  us  confidence  that  both  approaches 
are  essentially  correct. 


Summary 

In  spite  of  the  failure  of  ITT-GTC  to  deliver  their  devices,  we  made  substantial 
progress  in  analyzing  the  behavior  of  the  CDO  device  under  conditions  where  the 
source/drain  field  competes  with  the  gate  field  for  control  of  the  channel.  A  major 
accomplishment  was  the  creation  of  a  new  computer  code  which  accounts  for  both 
source/drain  and  gate  fields  simultaneously.  These  computer  simulations  reinforce 
our  belief  that  single-domain  operation  observed  in  the  experimental  devices  is  a 
result  of  loss  of  gate  control  of  the  channel. 

1 

However,  no  matter  how  sophisticated  the  computer  analysis,  it  can  never  replace  a 
working  experimental  device.  Unfortunately,  the  goal  of  experimental  verification 
was  not  achieved  in  this  project.  It  is  the  author’s  belief  that  the  contiguous  domain 
effect  is  real  and  can  be  observed  in  properly-designed  devices.  However,  the 
author  also  feels  that  he  has  personally  devoted  enough  time  and  effort  to  this  goal, 
and  that  any  further  development  work  should  be  carried  out  by  others.  He  is  eager 
to  share  his  thoughts  and  calculations  with  other  groups  who  might  be  inclined  to 
pick  up  this  reparch  in  the  future. 

The  author  wishes  to  express  his  sincere  appreciation  to  the  Air  Force  Office  of 
Scientific  Research,  and  in  particular  to  Dr.  Gerald  Witt,  for  their  support  of  this 
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in  this  project,  and  continued  to  support  this  work  at  a  low  level  after  the  original 
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Figure  7.  Electron  density  predicted  by  the  new  program  when  the  gate  is  placed 
in  close  proximity  to  the  channel  so  that  the  gate  field  dominates.  For  an  initial 
condition,  we  established  a  single  domain  in  the  channel  at  t-0.  but  by  10  ps  the 
channel  had  spontaneously  developed  multiple  domains,  and  the  pattern  had 
reached  steady  state  well  before  100  ps.  (Tlie  distribution  pictured  here  was 
observed  at  150  ps).  This  represents  the  CONTIGUOUS-DOMAIN  MODE. 
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I'igurc  8.  Frequency  of  oscillation  as  a  function  of  tlie  parameter  l/(ns  Xd  ), 
where  n,  is  the  average  electron  density  in  the  channel  and  Xd  is  the  separation 
between  the  centroid  of  channel  charge  and  the  gate  electrode.  The  cross, 
triangle,  star,  diamond,  and  circle  symbols  represents  results  from  the  original 
CDO  simulation  program  of  Fig.  4.  Tlic  large  black  circles  represent  results  from 
the  new  program  for  the  case  where  the  gate  field  dominates  over  the  source/drain 
field.  Agreement  is  very  good.  Tliis  gives  us  confidence  in  the  predictions  of 
BO'^1 1  programs. 
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grant  had  expired.  As  a  result  of  his  efforts,  two  graduate  students  were  able  to 
continue  their  research  and  obtain  their  Ph.D.  degrees  after  the  termination  of  the 
original  grant.  Both  arc  now  employed  in  the  U.S.  semiconductor  industry. 
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Outside  Interactions 

We  have  collaborated  with  the  ITT  Gallium  Arsenide  Technology  Center, 
Roanoke,  VA,  in  the  fabrication  and  characterization  of  the  initial  CDO  devices. 
The  devices  were  designed  and  masks  were  fabricated  at  Purdue.  ITT  conducted 
all  the  fabrication  operations  and  delivered  three  two-inch  CDO  wafers  to  Purdue 
on  23  February,  1988.  After  characterization,  we  determined  that  the  gate 
resistivity  of  the  original  devices  was  too  high,  and  ITT  agreed  to  fabricate  a 
second  lot  of  devices  for  delivery  in  the  spring  of  1991.  Due  to  problems  with  the 
processing  line  in  the  spring  and  summer  of  1991,  these  samples  were  never 
delivered. 

Over  the  years,  we  have  been  contacted  by  several  groups  who  have  had  some 
level  of  activity  in  investigating  this  device  concept.  These  groups  included  MIT, 
Jet  Propulsion  Laboratory,  Columbia  University,  and  Magnavox  (Torrance,  CA). 
Both  MIT  and  Columbia  stated  that  they  conducted  their  own  computer 
simulations,  and  both  reported  observing  microwave  oscillations  in  their 
simulations.  JPL  listed  a  research  program  on  this  d  ;vice  in  their  research 
summary  brochure  for  a  recent  year  (1990?),  but  this  program  is  not  now  believed 
to  be  active.  Magnavox  sent  a  representative  to  Purdue  several  years  ago  (1986?) 
to  discuss  fabricating  the  device,  but  we  subsequently  elected  to  collaborate  with 

nr. 

There  have  been  no  new  interactions  during  the  last  12  months. 
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CHAFEERl 

INTRODUCTION 


1.1  Introduction 

Since  solid  state  devices  were  introduced  into  microwave  area,  they  have  proven  to  be 
very  desirable  in  terms  of  reducing  the  circuit  size  and  cost,  enhancing  the  circuit 
performance,  and  increasing  the  reliability.  Consequently,  they  are  constantly  replacing 
microwave  tubes  and  continuing  to  open  up  new  applications  particularly  at  lower  power 
levels  [1].  In  many  cases,  they  are  the  heart  and  soul  of  numerous  microwave  systems. 
Especially  with  the  integrated  circuit  techniques  extending  to  the  fabrication  of  microwave 
circuits,  it  is  now  possible  to  reduce  greatly  the  size  of  microwave  circuits  so  that  complete 
microwave  systems  can  be  put  on  a  single  semiconductor  chip. 

The  microwave  power  generators  and  local  oscillators  are  essential  parts  of  the 
microwave  system.  Remarkable  successes  in  the  invention  and  improvement  of  the 
oscillators  have  been  achieved  over  the  last  30  years,  among  which  two  type*  of  oscillator 
arc  most  distinguishable,  namely  IMP  ATT  and  Gunn  diode  oscillators.  Various  practically 
used  oscillators  have  them  as  the  active  devices. 

To  meet  the  ever  growing  demands,  continuous  serrehes  for  new  devices  and  hew 
applications  are  still  going  on.  The  development  of  solid  state  nucrowave  sources  has  in 
general  been  of  two  kinds  [2].  First  of  all,  we  have  the  steady  progress  resulting  from 
pushing  the  upper  limiting  operation  frequency  of  the  existing  devices  and  improving  their 
compatibility  with  modem  VLSI  fabrication  technology  in  order  to  have  better  circuit 
performance,  low  cost  and  weight  Secondly,  we  have  the  more  erratic  progress  resulting 
from  the  discovery  of  new  devices  working  on  i^undamcntally  different  principles.  As  for 
the  second  category,  numerical  simulations  have  come  to  play  an  increasingly  important 
role  as  a  tool  for  exploring  new  concepts  and  designing  specific  devices  since  realistic 
modeling  of  nonlinear  behavior,  based  on  purely  analytical  techniques,  is  limited. 

The  current  research  project  is  on  a  potential  microwave  source  with  many  novel 
properties.  It  combines  the  transferred  electron  effect  in  n-type  GaAs  material  together  with 
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the  internal  two-dimensional  electrostatic  geometry  to  offer  high  frequency  operation  and 
high  frequency  tunability.  This  device  is  referred  to  as  the  Contiguous  Domain  Oscillator 
(CDO)  based  on  the  fact  that  a  contiguous  sequence  of  domains  can  be  generated  in  the 
channel.  The  frequency  is  determined  by  the  spacing  between  adjacent  domains  rather  than 
by  a  transit  time.  Therefore  high  frequency  operation  is  achievable  without  submicron 
technology.  Last  but  not  the  least,  it  has  a  planar  structure  and  can  be  easily  integrated  into 
the  microwave  subsystems  now  under  development. 


1.2  Comparative  Description  of  Conventional  Diodes  and  the  Contiguous  Domain 
Oscillator  (CDO) 

1.2. 1  Conventional  Diodes  ; 

By  conventional  diodes,  we  mean  IMPATT  and  Gunn  diodes.  The  IMP  ATT  diode  is  a 
representative  of  the  transit  time  device  family,  and  the  Gunn  diode  is  classified  as  a 
transferred-electron  device.  Although  the  charge  formation  mechanism  is  different  for  each 
type  of  device,  they  all  operate  using  the  effect  of  a  finite  carrier  transit  time  as  had  been 
done  with  microwave  electron-beam  devices  of  the  Klystron  and  traveling-wave  tube  type. 

Therefore  they  are  normally  operated  into  an  external  resonant  circuit  to  achieve  better 

! 

efficiency,  frequency  stability  and  noise. 

As  shown  in  Fig  1.1(a),  the  conventional  di^es  have  one-dimensional  geometry.  A 
high  electric  field,  which  is  above  threshold  for  the  formation  of  charge  packets,  is 
establisht  1  by  applying  a  voltage  Vq  across  the  ^o  terminals  (for  IMPATT  diode,  ^th  is 
the  maximum  electric  field  at  breakdown,  and  for  Gunn  diode,  ^ih  is  the  field  where  the 
field-velocity  curve  peaks).  Once  a  charge  packet  is  ge.aerated,  tne  field  due  to  the  space- 
charge  and  its  image  is  parallel  to  the  applied  field  which  significantiy  alters  the  total  electric 
field  outside  the  domain  to  a  level  below  |,j;  and  thereby  turns  off  the  domain  generation 
mechanism.  As  a  result,  only  one  domain  can  exist  in  the  channel  and  drift  across  the 
device.  By  the  time  the  charge  domain  is  removed  at  the  anode,  the  electric  field  is  back  up 
to  the  original  higher  value  to  cause  the  formation  of  a  new  charge  domain  at  the  cathode. 
And  the  process  of  domain  build-up,  propagation,  and  extinction  is  repeated.  In  this  way, 
transit-time  current  oscillations  are  produced  which  have  a  fundamental  frequency 

f,=  l/T  =  v,/L 


(l.l) 
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Figure  1.1  Schematic  representation  of  (a)  conventional  microwave  diodes,  and  (b)  the 
contiguous  domain  oscillator. 

In  conventional  devices,  internal  space  charge  alters  the  electrical  field 
throughout  the  device,  but  in  the  CDO  structure,  internal  space  charge  alten 
the  field  locally.  This  leads  to  the  fundamental  differences  in  operation  of 
the  two  classes  of  devices  [3]. 


where  L  is  the  drifting  length.  Since  the  domain  drifting  velocity  v  d  is  approximately  equal 
to  the  saturation  velocity,  the  operation  frequency  is  primarily  determined  by  L.  Once  the 
device  is  fabricated,  the  frequency  is  not  subject  to  significant  change. 

When  these  two- terminal  devices  are  connected  to  a  load  circuit,  Vq  is  no  longer 
constant  Their  one-dimensional  electrostatics  are  closely  coupled  to  the  voltage  excursions 
in  the  external  circuit  A  Umited  amount  of  electronic  frequency  tuning  is  often  achieved  by 
varying  circuit  reactance  using  a  varactor  diode. 

At  present,  all  conventional  diodes  are  produced  in  discrete  form  and  operate  in  an 
external  resonant  circuit.  This  technology  is  a  major  limitation  to  their  widespread  use  in 
compact,  lightweight,  low-cost  systems  that  require  a  higher  degree  of  integration. 
Moreover,  the  device  performance  depends  on  both  intrinsic  properties  prior  to  packaging 
and  its  impedance  matching  in  the  waveguide  circuit  through  the  packaging  parasitic  circuit 
elements.  Therefore,  the  application  of  monolithic  fabrication  techniques  to  the  oscillators  is 
desired.  However,  monolithic  fabrication  of  these  diodes  and  their  matching,  bias,  and 
other  passive  circuits  on  the  same  chip  for  a  fully  integrated  high  performance  microwave 
subsystem  is  not  a  easy  task.  Much  research  work  has  been  devoted  to  that  area  [4-7]. 

1 .2.2  Contiguous  Domain  Oscillator  (CDO) 

In  order  to  allow  a  contiguous  sequence  of  charge  domains  to  coexist  in  the  drifting 
channel,  it  i.s  necessary  to  modify  the  electrostatic  boundary  conditions  so  that  the  local 
field  of  eacn  domain  are  screened  from  the  rest  of  the  channel.  This  is  the  basic  difference 
between  the  conventional  oscillator  and  the  CDO. 

By  employing  a  resistive  gate  on  the  device,  it  is  possible  to  achieve  space  charge 
screening  provided  that  the  resistance  of  the  gate  is  in  the  right  range.  As  shown  in  Fig. 
1.1(b),  the  high  electric  field  provided  by  the  resistive  gate  voltage  Wqq  biases  the  device 
into  a  region  of  negative  differential  mobility  and  causes  charge  domain  to  form.  Tiie  image 
charge  of  the  domain  resides  on  the  resistive  gate,  so  that  nearly  all  of  the  space-charge 
field  is  normal  to  the  drift  direction.  Due  to  this  special  geometry,  the  presence  of  one 
charge  domain  does  not  alter  the  electric  field  in  the  remainder  of  the  device.  This  permits  a 
number  of  domains  to  coexist  in  the  drifting  channel  simultaneously. 

The  operation  of  the  CDO  is  shown  in  Fig  1.2.  The  barrier  between  the  resistive  gate 
and  the  channel  enforces  that  electrons  travel  along  the  channel  direction.  This  barrier  could 
be  a  physical  barrier  provided  by  a  high  band-gap  material  in  a  MOSFET  structure  or  a 
potential  barrier  provided  by  positive  biases  of  source/drain  with  respect  to  gatel/gate2  in  a 
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Figure  1.2  The  operation  of  the  CDO  device. 


MESFET  structure.  Since  the  QDO  has  four  tcrminaJs  and  three  applied  voltages,  it  gives 
more  control  over  the  formation,  transportation  and  extraction  of  charge  packets.  The 
resistive  gate  is  used  to  provide  a  constant  high  electric  field  Vgg  /  L  in  the  channel  and  to 
screen  the  space  charge.  The  function  of  the  source- to- gate  1  voltage  VsG  is  to  control  the 
electron  injection  level.  The  drain  junction  at  the  high  voltage  end  of  the  gate  is  biased  so  to 
collect  any  electrons  that  drift  down  the  channel.  With  high  enough  Vdg.  the  current 
coming  out  of  the  drain  is  insensitive  to  the  potential  on  the  drain,  making  the  device  look 
like  a  current  source.  Therefore,  no  external  resonant  circuit  is  required,  which  greatly 
simplifies  the  circuit  design. 

Fig  1.3(a)  shows  the  electron  density  under  the  condition  of  constant  insertion  rate 
(constant  Vsg)  as  indicated  by  the  original  transient  Time-Of-Flight  (TOF)  program  [S]. 
The  charge  is  distance  d  away  from  the  gate.  The  different  plots  are  at  different  times  in  the 
simulation  as  the  charge  front  drifts  from  the  source  end  of  the  channel  towards  the  drain. 
The  oscillation  stabilized  after  several  picoseconds.  Eventually,  the  oscillation  fills  up  the 
whole  channel,  continuously  traveling  from  the  source  to  the  drain.  Thus  the  oscillation  is 
governed  by  the  inherent  periodic  quasi-steady  state  transport  properties  of  the  carriers 
themselves. 

Frequency  tunability  is  realized  by  varying  source  to  gate  voltage  Vsg  which  controls 
charge  injection  into  the  channel.  When  the  electron  concentration  at  the  source  is  changed 
from  one  value  to  a  larger  value  for  a  period  of  time  and  then  remmed  to  the  original  value 
(Fig.  1.3(b)),the  oscillation  frequency  changes  accordingly.  This  is  a  strong  indication  that 
the  oscillation  frequency  can  be  electrically  tuned  instantaneously. 

Since  the  CDO  has  a  planar  structure,  its  integration  into  the  developing  monolithic 
microwave  system  is  straightforward. 


Figure  1.3(a)  The  time  evolution  of  electron  density  in  th'*  MOSh'ET  channel  as  electrons 
are  introduced  from  the  source,  llic  source  is  nxxlclcd  by  a  constant  surface 

electron  density  at  n,  eT.Ss  10’  ’  cm'^  from  y=0  to  y=5  jim.  d=500  A  and 
drift  field  vd  «  7000  V/cm  (8). 


Figure  1.3(b)  Effect  of  a  step  change  in  the  injection  level  on  the  oscillation 
(MOSFET).The  injection  level  was  initially  7.5x10”  but  was 
increased  to  1.5x10”  cm"2  after  100  ps  and  returned  to  7.5x10”  cm”  after 
200  ps  [8). 
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1 ,3  Tlicsis  Overview 

The  contiguous  dotnain  oscillator  has  great  promise  to  become  a  useful  microwave 
device  for  monolithic  millimeter-wave  integrated  circuit  (MMICs)  applications.  The 
ventlcation  of  its  cxi.stencc,  however,  is  not  a  simple  task.  Tliis  thesis  is  an  important  part 
of  the  continuing  effort  towards  the  understanding  and  realization  of  the  contiguous  domain 
oscillation  mode. 

This  research  project  has  gone  through  simulation,  design,  fabrication  and 
characterization  stages,  but  not  necessary  in  that  order.  Some  stages  have  been  gone 
through  several  times  for  further  investigation.  The  thesis  is  organized  in  such  a  way  that  it 
gives  reader  an  overall  view  of  the  project.  The  order  used  does  not  represent  the  actual 
time  sequence  of  the  project.  Some  previous  results  have  also  been  included  for 
completeness. 

Chapter  2  presents  the  steady-state  analysis  of  the  ideal  RG-FETs,  assuming  a  constant 
mobility.  The  analysis  provides  useful  information  about  the  regions  of  operation  in  the 
ideal  RG-FET  as  well  as  the  uniform  channel  condition  and  the  pinch-off  condition. 
Chapter  3  summarizes  the  transient  simulations  of  the  ideal  RG-FETs  with  the  full  velocity- 
field  curve.  Oscillations  ranging  from  20  GHz  to  60  GHz  are  observed.  The  dependence  of 
frequency  and  power  on  electron  density  in  the  channel  and  the  distance  between  gate  and 
channel  arc  presented.  Chapter  4  first  describes  the  fabrication  procedure  of  the  CDO 
devices  made  with  the  RG-MESFET  structure,  then  presents  the  FV  curve  and  microwave 
measurement  results.  Single  domain  oscillation  mode,  rather  than  contiguous  domain  mode 
was  observed  due  to  the  high  gate  resistance.  Chapter  4  also  presents  the  CrSiO  process 
development  for  the  purpose  of  modifying  the  single  domain  RG-MESFET  oscillator.  In 
order  to  further  understand  the  device,  a  simple  equivalent  circuit  model  is  established  in 
chapter  5  to  simulate  the  different  oscillation  modes  in  the  RG-PTET  devices.  Finally,  the 
conclusion  and  the  recommendation  for  the  future  work  are  given  in  Chapter  6. 
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CHAPTER  2 

STEADY  STATE  ANALYSIS  OF  THE  IDEAL  RG-FET 


2.1  Introduction 

In  considering  the  behavior  of  the  resistive-gate  MOSFET  device,  particularly  the 
CDO  application,  it  is  important  first  to  understand  the  carrier  behavior  inside  the  device 
in  the  DC  steady  state.  This  knowlege  provides  a  reference  frame  for  analyzing  the 
transient  bch.ivior,  and  can  indicate  the  proper  bias  conditions  for  contiguous  domain 
oscillation. 

When  the  gate-to-gatc  voltage  Vqg  set  to  zero  in  an  RG-FET,  it  becomes  a 
conventional  FET.  The  static  analysis  of  conventional  FETs  is  well  understood  [9].  With 
non-zero  Vco*  several  interesting  device  states  can  occur.  One  of  the  most  interesting  is 
the  uniform  channel  condition,  where  neither  the  channel  charge  per  unit  area  nor  the 
channel  electric  field  vary  from  source  to  drain.  Unlike  the  conventional  FET  device, 
where  the  uniform  channel  exists  only  when  the  drain-to-source  voltage  Vos  is  zero, 
Coen  and  Muller  have  shown  that  the  RG-MOSFET  can  have  a  uniform  channel  for  non¬ 
zero  VoSf  where  appreciable  drain  current  is  availablc[10,ll]. 

It  is  the  purpose  of  this  work  that  the  basic  Id  -  Vd  characteristics  of  the  RG-FET 
device  be  understood  as  well  as  the  free-carrier  distribution,  electric  field,  potential 
profile  along  the  channel  for  all  biases.  These  quantities  are  important  to  predict  the 
electric  characteristics  of  the  device  and  lead  to  an  easy  physical  interpretation  of  the 
modes  of  operation. 

Two  types  of  structure  arc  of  interest  for  the  simulation:  RG-MOSFET  structure  and 
RG-MESFET  structure.  Four  basic  assumptions  employed  for  the  following  work  are 
listed  here:  (1)  gradual  channel  approximation  [13]  is  applied,  which  states  that  the  rate  of 
change  of  the  electrostatic  variables  (potential,  electric  field,  etc.)  in  the  channel  direction 
are  relatively  smaller  than  the  rate  of  change  of  the  same  variables  in  the  vertical 
direction,  (2)  the  resistive  gate  is  a  perfect  linear  voltage  divider,  (3)  current  flow  is 
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parallel  to  the  channel  direction,  i.e.,  leakage  of  electrons  from  the  channel  to  the  gate  is 
negligible,  and  (4)  the  velocity  -  field  relation  is  linear.  Despite  the  above  assumptions, 
the  simulation  does  provide  useful  insight  into  the  basic  physical  behavior  of  the  device. 

To  model  the  RG-MOSFET,  we  extend  the  Pao-Sah  model  to  include  the  linear 
potential  drop  along  the  resistive  gate.  This  generalized  Pao-Sah  model  is  verified  by 
comparison  to  a  Time-Of-Flight  transient  simulation  for  the  RG-MOSFET.  We  also 
derive  analytic  expressions  for  the  Id-Vd  characteristics  which  are  sufficiently  accurate 
for  first-order  calculations. 

The  basic  I-V  characteristics  for  the  RG-MESFET  structure  are  also  examined.  This 
study  is  obtained  by  the  TOF  program.  The  uniform  channel  condition  and  the  pinch-off 
condition  are  derived  in  a  similar  way  as  for  the  RG-MOSFET. 


2.2  Generalized  Pao-Sah  Model  for  the  RG-MOSFET 

The  static  IV  characteristics  of  a  conventional  MOSFET  were  derived  by  Pao  and  Sah 

[12]  based  on  a  one-dimensional  model  subject  to  the  gradual  channel  approximation 

[13] .  The  same  approach  can  be  extended  to  compute  the  current- voltage  characteristics 
for  the  RG-MOSFET. 

Fig.  2.1  shows  a  cross-sectional  view  of  an  ideal  RG-MOSFET  structure,  indicating 
standard  terminals  and  dc  voltage  designations.  As  in  the  conventional  MOSFET,  carriers 
enter  the  structure  through  the  source  (S),  leave  through  the  drain  (D),  and  are  subject  to 
the  control  or  gating  action  of  the  gates  (G1  &  G2).  The  voltage  applied  to  the  source 
relative  to  the  neutral  substrate  is  Vs,  the  gate-1  voltage  relative  to  the  substrate  is  Vqi 
and  the  drain  voltage  relative  to  the  substrate  is  Vd.  The  voltage  applied  between  the 
gate-1  and  gate-2  is  Vgg-  Finally,  consistent  with  the  source  and  drain  designations,  the 
drain  current  Id  for  the  p-bulk  device  is  taken  to  be  positive  when  flowing  from  drain  to 
source  inside  the  device. 


2.2. 1  Basic  Equations 

Neglecting  the  contribution  from  holes,  the  drain  current  of  an  RG-MOSFET  can  be 
related  to  the  electron  quasi-Fermi  potential  gradient,  dvN  /dy,  and  the  sheet  charge 
density  Qn.  by: 


lD  =  -WQN(y)lin 


(2.1) 
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Figure  2.1  Cross-sectional  view  of  the  ideal  RG-MOSFET  structure  showing  the 
terminal  designations  and  standard  bias  condition. 


where  vn  =  (Ep  -  Fn  )  /  q,  W  is  the  channel  width,  and  iXn  is  the  electron  effective  surface 
mobility.  After  introducing  dimensionless  quantities  (see  Fig.  2.2)  and  integrating  with 
respect  to  y  on  both  sides  of  (2.1),  Id  can  be  rewritten  as: 


WMni^T 

L  q 


(2.2) 


Here, 


Uj(x)  =  [E(bulk)-E(x)]fkaT 


is  the  electrostatic  potential,  1 

Up  =  ( E.(bulk)  -  Ep  ]  /  kgT|=  In  ( N^/  n. ) 

l 

is  the  doping  parameter,  and 

C(y)  =  [Ep-Fj,]/kBT  =  v^yv.j. 


(2.3) 


(2.4) 


(2.5) 


is  the  quasi-Fermi-level  splitting,  all  normalized  to  the  normal  voltage  vt  =  ksT/q.  Uis  is 
the  Ui  value  at  oxide-semiconductor  interface  «x=0),  L  the  channel  length,  nj  the  intrinsic 
carrier  concentration,  Na  the  substrate  doping(assumed  uniform),  Us  the  source  voltage 
in  units  of  vt,  and  Ud  the  drain  voltage  in  units  of  vj. 
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Figure  2.2  Graphical  definition  of  dimensionless  variables  Uj  (x).  Up  and 


The  sheet  charge  density  is  expressed  as  [12] 

o  '  r  i 


pC(y))  I 


where 


V  ’  (2.7) 

is  the  intrinsic  Debye  length,  £$  is  the  semiconductor  dielectric  constant,  and 

F(UpUp,C(y))  =  V  c\u,-l )  +  e  "( (j.S) 

is  the  normalized  electric  field  obtained  fix>m  the  1-D  solution  of  Poisson's  equation. 

The  normalized  surface  potential  Uis  in  (2.6)  is  related  to  the  position  dependent  gate 

voltage  Vq  (y)  through  Kirchoff  s  voltage  law: 

k-T  e 

V  _ ril  J.  \ _ » _ TT/Tr  TT  1 


where 


Vo(y) = -2-  [  Ub  4-  sim(Ug)  nu^.  Up  0  j 

Vofy)  -  V^y)  -  V,,  -  v,,, + 


(2.10) 
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VpB  is  the  flat  band  voltage,  Sign(Uis)  the  sign  of  Uis,  and  Cqx  the  oxide  capacitance  per 
unit  area.  By  integrating  (2.2)  with  respect  to  an  arbitrary  the  normalized  distance 
along  the  channel  measured  from  the  source,  y/L,  may  be  related  to 


?(y) 

_ 

L“  Uq 

JQn(C)c1C 

Us 


(2.11) 


This  relation  is  not  valid  in  the  post-pinch-off  region. 


2.2.2  Algorithm 

Equations  (2.2),  (2.6),  (2.9),  (2.10)  and  (2.11)  constitute  a  complete  solution  for  the 
cunent-voltage  characteristics  of  the  RG-MOSFET,  but  they  must  be  solved  iteratively. 
Given  Voi,  Vgg.  and  Id,  the  drain  voltage  Vd  can  be  calculated  by  the  following  steps: 

(1)  At  the  source,  set  y=0  and  V'g  (0)=Vgi  -  Vfb- 

(2)  Guess  a  quasi-Fermi  potential  ^  at  that  y  position  and  solve  (2.9)  iteratively  for  the 
surface  potential  Uis- 

(3)  Substitute  the  calculated  Uis  into  (2.6)  and  numerically  integrate  for  the  inversion 
charge  Qn. 

(4)  Define  a  difference  function  A(y)  obtained  from  (2. 1 1)  as 

Wm  AT  cw 
A — JqnOcic 

u  u 

s  (2.12) 

and  calculate  A.  If  A  is  not  negligibly  small,  choose  a  different  ^  in  step  (2)  and 
repeat  steps  (2)  through  (4)  until  A  »  0. 

(5)  Qiange  position  y  and  calculate  V’g  (y)  using  (2.10).  Go  to  step  (2), 

When  Vgg  is  set  to  zero,  the  FV  characteristics  obtained  from  this  generalied  Pao-Sah 
(GPS)  model  are  identical  to  the  IV  characteristics  obtained  using  the  Pao-Sah  model. 
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2.2.3  Static  Simulation  Results 

All  the  results  shown  here  are  computed  for  an  n-channel  RG-MOSFET  device  with 
the  following  parameters:  Na  =  IxlO^^  cm  *3,  L  =2  ^m,  W  =  1  (im,  Pn  =  1500  cm2  /v-s, 
tox  =  500A,  ni  =  1.45x1010  cm-3,  gj  =  11.8  eq,  and  Eox  =  3.9  eq,  where  Eq  =  8.854x10-1'! 
F/cm. 


2.2.3. 1  rv  Curves  for  Different  Biases 

The  simulated  Id  -  Vd  characteristics  using  the  generalized  Pao-Sah  model  are 
presented  in  Fig.  2.3  as  solid  lines.  Fig.  2.3(a)  is  the  IV  curves  for  fixed  Vs=lV,  Vgi=6V 
and  different  gate-to-gate  voltage  Vgg  frotn  0  to  8V,  corresponding  to  an  average  electric 
field  from  0  to  40,000V/cm.  One  can  see  that  the  Id  -  Vd  curves  of  the  RG-MOSFET 
have  the  same  general  shape  as  the  conventional  MOSFET  (Vgg=0).  except  that  the 
pinch-off  voltage  Vdp  and  the  saturation  current  iDsat  are  increased  due  to  the  application 
of  gate-to-gate  voltage  Vgg-  We  can  gain  insight  into  the  physical  processes  taking  place 
in  the  device  by  examining  the  condition  along  the  channel  for  one  specific  Vgg  value. 
This  is  done  in  Fig.  2.4  for  Vgg  =  5V.  The  solution  for  Vco  *  0  are  shown  in  Fig.  2.5, 
for  comparison.  From  Fig.  2.4,  we  first  note  that  at  low  values  of  drain  current,  a  nonzero 
gate-to-gate  voltage  Vgg  causes  the  inversion  charge  density  to  increase  along  the 
channel.  Therefore,  the  device  has  a  nonunifonn  channel  at  Vds=0.  If  Vg  ,  Vgi,  and  Vgg 
are  fixed  (Fig.  2.4(a)),  the  drain  current  increases  with  increasing  drain  bias  Vd, 
corresponding  to  tlie  line  from  point  A  to  point  U.  We  shall  call  this  the  "pre-uniform" 
region,  where  Vd  is  less  than  the  drain  voltage  Vdu  corresponding  to  the  uniform 
channel  condition.  In  the  pre-uniform  region,  the  diffusion  current  is  determined  by  the 
slope  of  the  electron  density,  which  is  in  the  opposite  direction  to  the  drift  current,  as 
shown  in  Fig.  2.6  (dots).  The  electron  density  at  the  drain  is  greater  than  at  the  source  and 
is  reduced  by  the  applied  drain  voltage,  as  shown  in  Fig.  2.4(b).  Fig.  2.4(c)  shows  that  the 
electric  field  is  nearly  uniform  along  the  channel. 


Figure  2.3  Id*Vd  characterisdcs  for  an  RG-MOSFET  given  by  the  generalized  Pao- 
Sah  model  (lines).  Bias  conditions  are:  (a)  Vs  =  IV  and  Vgi  =  6V,  (b)  Vs 
=  IV,  Vgg  =  5V.  Stars  indicate  pinch-off  voltages  calculated  from  (2.21), 
triangles  indicate  the  uniform  channel  condition  calculated  from  (2.20), 
and  circles  are  currents  calculated  by  a  one-dimensional  TOF  transient 
simulation[6]. 
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When  the  electron  density  at  the  source  is  juit  equal  to  the  electron  density  at  the 
drain,  the  RG-MOSFET  has  a  uniform  channel  (point  U).  The  electric  field  in  the  channel 
is  now  uniform,  and  the  diffusion  current  is  zero  (Fig.  2.6,  open  diamonds).  Notice  that  in 
the  vicinity  of  Vdu.  the  surface  channel  acts  more  or  less  like  a  simple  resistor. 

For  Vd>Vdu.  the  device  is  in  the  "post-uniform"  region.  The  electron  density  at  the 
drain  decreases,  reaching  zero  when  Vd=Vdp.  This  corresponds  to  the  line  from  point  U 
to  point  P  in  Fig.  2.4(a).  The  electron  density  at  the  drain  decreases,  reaching  zero  when 
Vd=Vdp.  This  is  referred  to  as  the  "pinch-off  condition.  The  region  where  Vp  >Vdp  is 
designated  as  "pinch-off  region.  In  this  region,  the  electric  field  at  the  source  is  pinned  at 
its  maximum.  The  diffusion  current  is  now  positive  and  increases  dramatically  at  the 
drain  ( Fig.  2.6,  open  circles). 

The  reduction  of  the  channel  charge  caused  by  increasing  the  drain  bias  is 
counteracted  by  the  increase  in  channel  charge  caused  by  the  linear  increase  of  the  gate 
bias.  For  this  reason,  pinch-off  at  the  drain  is  delayed  and  the  saturated  drain  current  is 
considerably  larger  than  in  the  MOSFET.  The  bigger  the  Vgg«  the  larger  the  drain  current 
Id  at  the  same  drain  bias. 

Comparing  Fig.  2.4(c)  and  Fig.  2.5(c),  one  can  see  that  the  nonzero  Vgg  raises  the 
overall  electric  field  in  the  channel.  Therefore,  high  electric  field  disfribution  across  the 
channel  can  be  achieved  in  the  RG-MOSFET  structure. 

The  inverted  region  of  the  conventional  MOSFET  is  similar  to  the  pre-uniform  region 
of  the  RG-MOSEET  if  electron  density  distribution  is  examined.  In  both  case,  the  drain 
has  more  electrons  than  the  source  does.  However,  the  drain  current  is  negative  in 
conventional  MOSFET. 

Fig.  2.3(b)  are  the  IV  curves  for  fixed  Vs=lV  and  Vgg=5V,  and  different  gatel 
voltage  Vgi  from  4V  to  IIV.  As  Vgi  increases,  the  channel  electron  density  increases, 
and  drain  current  increases  proportionally. 

In  Fig.  2.3(a).  the  results  of  a  one-dimensional  finite  difference  simulation[6]  for  Vgg 
=  5V  are  shown  as  circles.  The  finite  difference  simulation  not  only  predicts  the  same 
steady-state  current  as  the  generalized  Pao-Sah  model,  but  gives  the  same  electron 
density,  electric  field,  and  channel  potential  distributions  as  well. 
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Figure  2.4  (a)  Id-Vd  curve,  (b)  surface  electron  density  nj,  (c)  electric  field  ^y,  and 

(d)  Fermi  potential  vn  along  the  channel  for  the  RG-MOSFET  of  Fig.  2.3, 
biased  at  Vs  =  IV,  Vqi  =  6V,  and  Vgg  =  5V,  calculated  using  the 
generalized  Pao-Sah  model.  Dots  in  (a)  correspond  to  the  lines  in  (b),  (c) 
and  (d).  "U"  and  "P"  mark  the  uniform  channel  and  pinch-off  conditions, 
respectively. 
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igure  2.5 
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(a)  Id-Vd  curve,  (b)  surface  electron  density  nj,  (c)  electric  field  ^y,  and 
(d)  Fermi  potential  vn  along  the  channel  for  a  conventional  MOSFET 
(identical  to  the  RG-MOSFET  of  Fig.  2.3  except  with  Vgg  -  0),  biased  at 
Vs  =  IV  and  Vgi  =  6V,  calculated  using  the  standard  Pao-Sah  model. 
Dots  in  (a)  coirespond  to  the  lines  in  (b),  (c)  and  (d). 
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Figure  2.6  Diffusion  current  components  for  the  RG-MOSFET  biased  at  Vs=lV, 
Vgi=6V  and  Vgg=5V.  A,  U,  and  P  correspond  to  Fig.  2.4(a). 


2.2.3.2  Electric  Field  Distribution  in  the  RG-MOSFET 

One  has  to  remember  that  the  generalized  Pao-Sah  simulation  is  based  on  the  gradual 
channel  approximation.  The  simulation  is  only  valid  when  »  ^y-  Fig-  2.7  compares 
the  channel  electrical  field  distribution  ^y(y)  with  the  vertical  electrical  field  distribution 
in  the  saturation  region  for  different  values  of  Vgi  and  Vgg- 

Fig.  2.7  shows  that  the  gradual  channel  approximation  (GCA)  is  violated  at  the  drain 
since  the  channel  is  first  pinched-off  there  when  lD=lDsat-  This  violation  also  occurs  for 
the  conventional  MOSFET.  However,  it  is  more  series  for  non-zero  Vqg-  For  fixed  gate- 
to-gate  voltage,  the  violation  to  the  GCA  occurs  when  Voi  is  small  (Fig.  2.7(b)). 
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Effect  of  (a)  gate-to-gate  bias  Vgg,  and  (b)  gatel  bias  Vgi  on  the 
distribution  of  electrical  fields  along  the  channel  (^y)  and  perpendicular  to 
the  channel  (£x).  Id  “Idshl 
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2.2.4  Simple  Analytic  Mtxicl  Tor  tl.c  RG-MOSFET 


2.2.4. 1  Unifonn  Channel  Condition  and  Pinch  off  Condition 

llie  analytical  equations  used  to  desenbe  the  conventional  MOSFET  can  be  applied  to 
the  resistive  gate  device,  provided  that  they  arc  modified  to  account  for  the  potential  diop 
along  the  gate: 

W  =  +  (2,13) 

In  strong  inversion,  the  surface  potential  is  2<)F'*^''N(y),  where  vNfy)  is  the  quasi  Fermi 
level  splitting  and  iftp  *  koT/q  ln(NA/nj)  is  the  Fermi  potential  of  the  substrate.  Therefore, 
the  sheet  electron  density  is; 

"tly)  =  -q-  [  y  ^  (2  14) 


where 


(2.15) 

is  the  body  factor.  Expanding  the  square  root  term  of  (2.14)  around  vn=Vs  [15],  we  have: 


_  _ _  V  (y)  •  V 

/^y)  +  2(5,.  =  /Vj  +2(5,  + 

2yVj.  +  2.5, 

Equation  (2.14)  becomes 

".(y)  -  %  1  V„,-  [-  .v^ .  (1  *y)  (v„(y) .  V,)  J 


where 


V~v^ 


(2.16) 


(2.17) 


(2.18) 


Vt  =  2vv,>/2v;/?^:.v3>v,„ 


T  's  y-  ov^F  ’s  tn  (2.19) 

is  the  threshold  voltage,  defined  as  the  gatc-1  voltage  relative  to  the  bulk  which  cau.sc.s  a 
depletion  inversion  transition  at  the  source  end  of  the  channel. 

A  uniform  channel  is  obtained  when  the  sheet  electron  density  and  the  electric  field  ,ire 
independent  of  y.  If  v^ly)  is  taken  to  lx;  a  nearly  linear  function  of  position,  v-sity)  » 

(y/l-).  the  uniform  channel  condition  is  derived  from  (2. 17)  { lO.1 1 ): 

V..., 


V  B  V  4-  '**, 

18/  $  I^Y 


(2.2(9 
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Vdu  values  calculated  using  (2.20)  are  indicated  by  solid  triangles  in  Fig.  2.3.  They  agree 
very  well  with  simulated  values.  Since  the  above  derivation  docs  not  involve  the 
characteristics  of  the  velocity-field  curve,  this  uniform  channel  condition  is  universally 
valid. 

If  the  drain  current  saturation  is  caused  by  the  pinch-off  at  the  drain,  Vdp  is  found  by 
setting  y=L  and  ns(L)=0  in  (2.17), 


'DP 


V(L) 


V  +  V  -  V 

^01  ^  ^GO 

1+7 


(2.21) 


Equation  (2.21)  has  a  ver  /  clear  physical  meaning.  Recall  that  in  the  post-uniform  region, 
the  RG-MOSFET  acts  more  or  less  like  a  MOSFET  with  a  pinch-off  voltage  equal  to 
Vs+(Vgi-Vt)/(1+Y).  But  before  this  region,  Vds  =  Vgg  /(1+Y)  has  to  be  applied  to 
achieve  the  uniform  channel.  Therefore,  Vdp  of  the  RG-MOSFET  is  bigger  than  Vdp  of 
the  MOSFET  by  an  amount  equal  to  Vgg/(1+Y)-  (2.21)  is  used  to  calculate  the  pinch-off 
voltages  in  Fig.  2.3;  these  values  arc  marked  by  stars  in  the  figure.  The  calculated  pinch-off 
voltage  agree  well  with  the  onset  of  current  saturation  given  by  the  GPS  model. 


2.2.4.2  Regions  of  Operation  of  the  RG-MOSFET 

As  discussed  before,  the  RG-MOSFET  has  three  regions  of  operation  as  summarized  in 
Fig.  2.8.  For  simplicity,  yis  set  to  zero.  The  U  line  indicates  a  uniform  channel  condition: 
Vds  =Vgg.  The  P  line  indicates  a  pinch-off  condition:  Vds=Vgi  -Vy  +Vgg-  All  devices 
experience  these  three  regions  when  the  drain  voltage  increases.  When  Vgg  is  fixed  (Fig. 
2.8(b)),  non-zero  drain  current  only  occun  for  Vgi>Vt.  When  Vgi  is  fixed  at  any  value 
greater  than  Vy  (Fig.  2.8(c)),  the  U  line  and  P  line  arc  parallel  to  each  other.  The  pre- 
uniform  region  docs  not  exist  when  VGG=0t  and  the  uniform  channel  occurs  at  Vds=0. 

Since  the  uniform  channel  condition  is  univcnally  valid  for  all  shapes  of  velocity-field 
curves,  the  U  line  is  important  to  the  identification  of  the  microwave  oscillation  in  practical 
RG-MOSFET  devices,  provided  that  Vdu/L  is  greater  than  the  threshold  field  for  negative 
differential  mobility. 


j 
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Figure  2.8 


^faG 


constant 


(b) 


(c) 


Mode  diagram  illustrating  the  regions  of  operations  (y=0).  I  is  the  Pre¬ 
uniform  region,  11  is  the  Post-uniform  region,  and  III  is  the  pinch-off 
region. 

(a)  Coordinates  used,  (b)  Vgg  is  fixed, (c)  Vci  is  fixed. 

U  Une:  Vds(U)=Vgg 
P  line:  Vds(P)=Vgi-Vim-Vgg 
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2.2.4. 3  Simplified  Analytic  Model 

An  analytic  expression  for  the  drain  current  can  be  obtained  by  neglecting  diffusion.  By 
analogy  with  the  conventional  MOSFET,  we  may  write: 

Wu_  k  dv„  ^ 


•d— 

0  (2.22) 

where  VN  is  the  channel  Fermi  potential.  From  (2.17),  the  sheet  charge  density  is 

Q„(y)  =  •  Co.  [  Vo(y)  -  -  (l+D  ( v„(y)  -  Vj)  ] 

=  ■  C..  [  Vo,  -  -  (14D  ( v„(y)  -  V,)  *  Voof  ] 

Substituting  (2,23)  into  (2.22)  and  integrating  the  right  hand  side  from  the  source  to  the 
drain: 

I  +  f  .X®.  y  dy  ] 

L  oi  2  DS  J  L  ^  dy 

0  (2.24) 

The  first  two  terms  corresponding  to  the  familar  square  law  equation  for  the  MOSFET 

[34].  The  last  term  can  be  evaluated  using  the  chain  rule. 


=  ^[LVjj-Jv^(y)dy] 

0  (2.25) 

As  shown  in  Fig.  2.4(d),  for  Vd  ^  Vdu  the  channel  Fermi  potential  VN(y)  can  be 
approximated  as  a  linear  function  of  position: 


®  ^DS  ^  ^  ^DU 

For  Vd  >  Vdu.  we  approximate  vN(y)  with  a  second-order  polynomial: 


(2.26a) 


v„(y)-v,+Vo„(i)*(v„.Vo„)(i)*  .Vo>Vbu 

Figure  2.9  shows  a  comparison  between  the  actual  potential  given  by  the  GPS  model  and 
the  approximate  potentials  given  by  (2.26),  Agreement  is  satisfactory,  except  at  drain 
biases  very  near  pinch-off,  Vd  *  Vdp. 

Using  (2.26),  (2.25)  can  be  expressed  as 
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Figure  2.9  Channel  Fermi  potential  given  by  the  GPS  model  (circles)  and  by  the 
analytical  equations  Gines)  for  the  device  of  Fig.  2.3,  biased  at  Vs  =  IV, 
Vqi  =  6V,  and  Vqg  =  5V.  The  analytical  equations  (2.26a&b)  provide 
good  agreement  with  the  GPS  model  except  when  the  drain  bias  is  close  to 
pinch-off. 


-  •! 


JLv  V 

2  GG  ^DS 


V  S  V 
*  D  DU 


0  ^  I  2.y  Y  . J.V 

3  '^GO  '^DS  6  CG 

Therefore,  the  drain  current  expression  (2.24)  becomes 

f  iy  V 

2  CO  '^DS 

,RG  _  ^MOS  ^\lnCox  J 

■  L  '9  1 

^  ^  \r  \r  I  \r 


—  V  V  --^V  V  V  >v 

a  '^GO '^DS  fi  ''gG  '^DU 


(2.27) 


J-V  V 

2  ^GG  '^DS 

V  ^  V 
’  DU 

■3  ^GG^D3  " 

Xv  V 

6  '^GG  ^DU 

V  >  V 
*  ’^D  DU 

(2.28) 

-Vt)Vds- 

1 - 1 

> 

(2.29) 

where 


Equation  (2.28)  reduces  to  the  familiar  MOSFET  current  equation  (2.29)  when  Vgg  is 
set  to  zero.  The  drain  current  given  by  (2.28)  reaches  a  maximum  when  V©  equals  Vosat 
given  by 
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V  +— V  -V 
V  -X/  ,  ^01^3^00 

°***  ^  (2.30) 

Vosat  is  slightly  smaller  than  the  pinch-off  voltage  given  by  (2.21).  For  Vd  >  Vcsat.  the 
drain  current  is 


,80  _  WivCox  <''01*3''00‘''t>  1.,  „  , 

‘O'"  ■  2L  L  (1+7)  ■3''®''ouJ  P3,J 

Figure  2.10  shows  a  comparison  between  the  current  predicted  by  the  analytical  model  of 
(2.28)  -  (2.31)  and  the  generalized  Pao-Sah  model.  Agreement  is  quite  reasonable  over  the 
whole  range  of  bias  voltages  investigated,  indicating  that  the  analytical  model  is  adequate 
for  first-order  calculations. 

Note  that  (2.29)  could  be  expanded  to  the  standard  "bulk  charge"  MOSFET  -equation 
[34]  if  (2.14)  were  used  instead  of  (2.17)  in  deriving  (2.22): 

^os_  . 


(2<l.pW5) 


(2.32) 


where 


-  2<|>f  +  VpB  +  Vj  (2.33) 

The  drain  current  for  a  RG-MOSFET  is  still  in  the  same  form  as  in  Equation  (2.28). 
With  this  "bulk  charge"  model,  reaches  a  maximum  when  V©  equals  Vosai  given  by: 

VDM.-Vs+[''oi-V™+fv<^+V,.^vJ*2V„CVo,-V„*fv<^] 
the  corresponding  uniform  channel  condition  is: 

Vw,  -  V  Voo-  v,+  V„-7'^+2V„(Voo-'',+  2‘I’f  +  ''s)  n  351 


where 


Vr-v'2V,(2<Dp+V5) 


(2.35) 


(2.36) 


and  the  pinch-off  voltage  is: 


VnP  =  V,-2<D,-H 


■  +  V„,+  V. 


cr  ’CC  ’FB 


“  '^S  '*■  ^gg'*'  -v/^  +  ^'^o^gi'^fb'^  ^gg\ 


5  L  Gl  FB  GG  OV  O  O'  GI  FB  GG'J  ^2  37) 

Fig.  2.11  compares  Vdp  calculated  from  Eqn  (2.22)  and  (2.37)  for  different  doping 
densities.  The  approximation  of  (2.16)  is  fairly  good  for  the  doping  densities  and  gate-to- 
gate  voltages  investigated. 


Drain  Cuirent  (mA)  Drain  Current  (mA) 
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Figure  2.10  Id*Vd  curves  given  by  the  GPS  model  (circles)  and  the  analytic  model 
(lines)  for  the  device  of  Fig.  2.3.  Bias  conditions  arc:  (a)  Vs  =  1 V  and  Vqi 
=  6V,  (b)  Vs  a  IV,  Vgg  =  5V.  Agreement  is  good,  indicating  the  analytical 
model  can  b;  safely  used  for  first-order  calculations. 
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Figure  2.11  Comparison  of  Vdp  calculations  with  and  without  the  approximation  of 
(2.16).  Lines  are  calculated  from  (2.21),  with  the  approximation,  x's  are 
calculated  from  (2.37),  without  the  approximation.  Vs=lV,  Vgi=6V,  and 
Vi=1.56V. 


2.3  Simulation  Results  for  the  RG-MESFET 

Although  different  in  structure  and  threshold  voltage,  the  operation  of  the  RG*MESFET 
is  similar  to  the  RG-MOSFET.  The  conventional  coordinate  system  is  used  as  shown  in 


Figure  2. 12  Cross  section  view  of  the  ideal  RG-MESFET  structure  showing  the  terminal 
designations  and  standard  bias  condition. 
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Fig.  2.12.  The  RG-MESFET  has  a  doped  channel  with  length  L,  width  W  and  a  doping 
density  No- 

2.3. 1 IV  Curves  for  the  RG-MESFET 

All  the  results  shown  here  are  for  an  n-channel  RG-MESFET  device  with  the  following 
parameters:  Nd=1.2  x  10^"^  cm*3,  L=2pim,  W=ljim,  pin  =5200  cm^  /V-s,  Tepi  =1500A, 
and  Es  =  12.9  £q.  The  simulations  are  done  using  the  Time-of-Flight  program  (for  a  detailed 
description  of  the  TOF  program,  see  Chapter  3). 

The  simulated  Id  -  Vd  characteristics  for  Vs  =0.5 V,  Vqi  =  0  and  different  gate-to-gate 
bias  Vgg  arc  presented  in  Fig.  2.13(a).  The  curves  of  non- zero  Vqg  are  similar  in  shape  to 
the  conventional  MESFET  (Vgg=0).  only  the  uniform  channel  and  pinch-off  voltages  are 
increased  due  to  the  same  reason  as  discussed  for  RG-MOSFET  in  section  2.2.3. 1.  The 
current  does  not  go  to  zero  for  non-zero  Vgg  due  to  the  assumption  in  the  program  that 
restrict  the  current  flow  only  in  the  channel  direction.  In  the  region  where  Vd<2V,  the 
Schottky  diode  is  forward  biased,  causing  a  large  amount  of  electrons  to  flow  from  the 
channel  to  the  gate.  The  simulation  is  not  valid  anymore.  When  Vd  <  Vs,  the  drain  current 
is  negative  for  Vgg=0. 

With  Vgg  fixed  at  2V,  the  IV  curves  for  Vgi=0  and  different  Vs  are  shown  in  2. 13(b). 
Increasing  Vs  reduces  the  number  of  electrons  at  the  source.  Consequently,  the  average 
electron  density  in  the  channel  is  decreasing. 

Fig.  2.14  and  2.15  show  the  comparison  between  the  RG-MESFET  and  the 
conventional  MESFET.  In  the  conventional  MESFET  (Fig.  2.14(b)),  uniform  channel 
occurs  only  when  Vd=Vs,  and  the  uniform  channel  current  is  zero.  In  the  case  of  Vd  < 
Vs,  the  drain  has  more  electrons  than  the  source,  and  electrons  flow  from  the  drain  to  the 
source,  resulting  in  negative  drain  current.  With  non-zero  Vgg  (Fig-  2.14(a)),  however.  Id 
is  positive  even  when  the  drain  voltage  is  smaller  than  the  source  voltage.  In  other  words, 
the  current  at  Vd=Vs  for  RG-MESFET  is  not  zero,  as  pointed  by  the  U  point  in  Fig. 
2.14(a). 

Again,  the  non  zero  Vgg  raises  the  overall  electric  field  in  the  channel  (Fig.  2.14(c)  & 
Fig.  2.15(c)).  High  electric  field  distribution  is  achievable  with  the  RG-MESFET  structure. 

So,  the  rv  curves,  electron  density,  electric  field  and  potential  distributions  for  an  ideal 
RG-MESFET  are  similar  to  that  of  the  ideal  RG-MOSFET.  There  are  also  three  regions  of 
operation  divided  by  the  uniform  channel  condition  and  the  pinch-off  condition.  However, 
unlike  the  RG-MOSFET,  where  the  maximum  electron  density  at  specific  point  in  the 


I 


channel  is  determined  by  the  gate  voltage  at  that  point,  the  maximum  electron  density  in  an 
RG-MESFET  channel  is  limited  by  the  product  of  channel  doping  density  and  the  channel 
thiclmess  NoTepi. 
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Figure  2.13  In-Vo  characteristics  for  an  RG-MESFET  given  by  the  TOF  progra.Ti 
(lines).Bias  conditions  are:  (a)  Vs  =  0.5V  and  Vqi  =  OV,  (b)  Vqi  =  OV, 
Vgg  =  2V.  Stars  indicate  pinch-off  voltages  calculated  from  (2.47), 
triangles  indicate  the  uniform  channel  condition  calculated  from  (2.46). 
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Figure  2.15 
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2.3.2  Uniform  Channel  Condition  and  Pinch-off  Condition  for  the  RG-MESFET 

Assuming  that  the  undepleted  channel  is  neutral,  the  position  dependent  sheet  electron 
density  in  a  doped  RG-MESFET  is: 

ns(y)  =  Nq  [x^(y)  -  Xj(y)] 

where  x<i  is  the  depletion  width  of  the  Schottky  junction,  and  Tepi  -  Xd2  is  the  depiction 
width  in  the  channel  due  to  channel-substrate  junction  (Fig.  2.16(a)).  If  the  substrate  is 
semi-insulating, 


^d2  “  ^epi 


(2.39) 


Xd  can  be  calculated  by  the  abrupt  p-H  n  junction  expression  [16]: 

y  (2.40) 

j  where  VN(y)  is  the  channel  Fermi  potential,  relative  to  the  Fermi  level  in  the  substrate,  and 

!  E„(channel)  -  F^, 

V,..  =  - - — 2^ 

bi  '*'Bn  q  (2.41) 

is  the  built-in  potential  for  a  Schottky  junction,  and  4)Bn  is  the  barrier  height  (Fig.  2.16(b)). 


Figure  2.16 


Cross  sectional  view  and  simple  band  diagram  for  a  MESFET  structure. 


For  uniform  ns,  the  quantity  inside  the  square  root  has  to  be  constant,  i.e.. 


VN(y)-VG(y)  =  const  ^^^2) 

When  the  uniform  channel  forms,  the  channel  Fermi  potential  is  a  luiear  function  of  the 
position: 


Vj^(y)  =  Vj+  Vjjg  — 

Since 

VoCy)  =  Vqj  +  Vqq  ^ 

the  uniform  channel  condition  is: 


i,e.. 


V  =  V 


V  =V  +  V 
''du  ''s  ’'gg 


(2.43) 


(2.44) 


(2.45) 


(2.46) 

Equation  (2.46)  is  similar  to  the  uniform  channel  condition  for  an  RG-MOSFET  (Equation 
(2.20)).  Values  of  Vdu  calculated  from  (2.46)  are  indicated  by  the  triangles  in  Fig.  2.13. 
The  pinch-off  voltage  can  be  obtained  by  setting  y=L  and  xj  =Tepi  in  (2.40): 


V  —V  a.  V  I  ^ 

^DP~  ''Gr  epi'  bi 

Vdp  values  calculated  from  (2.47)  are  shown  as  stars  in  Fig.  2.13. 


(2.47) 


2.4  Gjnclusion 

By  taking  into  account  the  gatc-to-gate  bias  Vgg.  the  steady-state  simulations  of  the 
long  channel  RG-FET  devices  with  constant  mobility  have  been  performed.  The  simulation 
from  both  generalized  Pao-Sah  model  and  Time-of-Flight  program  indicated  that  the  Id-Vq 
curve  of  the  RG-FET  have  the  same  general  shape  as  that  of  conventional  FET's. 
However,  internally,  the  operation  of  the  RG-MOSFET  experiences  three  regions  for 
positive  current.  They  are  divided  by  the  uniform  channel  condition  and  the  pinch-off 
condition.  The  pre-uniform  region  is  solely  caused  by  the  non-zero  Vgg-  The  post-uniform 
region  and  the  pinch-off  regions  are  very  simil:.r  to  the  conventional  FET’s  ohmic  and 
pinch-off  regions. 


^5 


With  non  zero  Vf.f;,  hif^hirr  cicctnc  field  values  arc  obtained  across  the  channel, 
a  '.uinin)j  constant  inobiliiy.  'llic  CIX)  operation  requires  high  electric  field  in  the  channel 
With  iraliMic  schx  ity  field  curve  of  GaAs.  charge  oscillation  is  p<n  ubic  once  the  electric 
field  djstributK'n  exceeds  the  tlursholii  field  Ijij,  in  a  resistive  gate,  MiT. 

I'he  unifonn  channel  condition  ts  universally  true  in  spite  of  the  shape  of  the  vchxity- 
ficld  curve.  This  condition  can  lx-  u'-d  as  an  initial  biasing  point  for  obtaining  the 
contiguous  donaain  oscillation. 
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CHAPTERS 

TRANSIENT  ANALYSIS  OF  THE  IDEAL  FET  DEVICE 


This  chapter  summarizes  the  large-signal  operation  characteristics  of  the  ideal 
contiguous  domain  oscillator,  as  obtained  by  the  original  Time-Of-Flight  (TOF)  program. 
This  TOF  program  was  developed  by  J.  A.  Cooper  in  1980  and  was  used  to  simulate  the 
high-field  velocity  of  electrons  in  inversion  layers  [14,  17].  The  detailed  procedure  of  the 
TOF  program  is  described  in  section  3.1.  Section  3.2  and  3.3  summarize  the  TOF 
simulation  results  for  both  RG-MOSFET  and  RG-MESFET  contiguous  domain  oscillators, 
where  the  effect  of  gate  control  and  the  sheet  electron  density  in  the  channel  on  the 
oscillation  frequencies  and  amplitudes  are  presented. 


3. 1  Description  of  the  Original  TOF  Program 

3.1.1  Basic  Assumptions 

The  basic  assumptions  used  in  the  TOF  program  are:  (1)  ideal  gate  (perfect  linear 
voltage  divider),  and  (2)  one-dimensional  current  flow  along  the  channel.  The  full 
velocity-field  curve  v  -  ^  is  used  in  the  simulation  in  order  to  observe  the  oscillation 
behavior. 

Although  different  in  structure,  both  the  RG-MOSFET  and  RG-MESFET  can  be 
simulated  with  the  TOF  program,  which  uses  a  finite  difference  technique  to  solve  the 
continuity  equation  for  the  electron  concentration  per  unit  width  ni(i)  along  the  channel. 
Fig.  3.1  shows  the  cross  section,  charge  density,  and  energy  band  diagram  along  th’ 
direction  perpendicular  to  the  channel  for  both  MOSFET  and  MESFET. 

It  can  be  seen  that  in  MOSFET  structure,  the  electrons  (minority  carriers)  are  confined 
to  the  interface  within  a  very  narrow  layer.  This  allows  the  removal  of  the  dimension  into 
the  substrate.  Surf.ice  electron  density  n,  is  easily  defined  and  so  is  the  electron  density  per 
unit  width:  ni«n,Ay.  In  the  MESFET  structure,  electrons  (majority  carriers)  move  along  the 
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undepleted  portion  of  the  channel  (N)  region.  The  assumptions  in  the  TOF  program  for  the 
MESFET  arc  that  the  electrons  occupy  the  undepleted  portion  of  the  epi-layer  with  a 
volume  density  equal  to  the  doping  density.  This  way,  the  electron  density  per  unit  width 
can  be  defined  as  ni  =  ND(xd2-Xd)Ay  for  each  grid  point  positioned  at  the  center  of  the 
charge  layer  Xm(i).  In  reality,  the  doped  channel  formed  by  ion-implantation  exhibits  a 
Gaussian  distribution.  One  can  closely  approximate  ion  implantation  profile  by  using  an 
effective  uniform  profile  [35]. 


Figure  3.1  (a)  Cross  section,  (b)  charge  density,  and  (c)  potential  from  gate  to  substrate 

for  MOSFET  and  MESFET  structures. 
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Having  made  the  above  assumptions,  the  numerical  treatments  of  MOSFET  and 
MESFET  really  are  the  same,  in  spite  of  the  fact  that  the  electrostatics  are  different  for  the 
two  structures.  Fig.  3.2  shows  the  equivalent  circuit  of  the  simulated  device  coiresponding 
to  the  cross  section  in  Fig.  2.1  (MOSFET)  and  Fig.  2.12  (MESFET).  The  gate  is 
represented  by  a  series  of  ideal  voltage  source.  The  channel  is  modeled  by  a  series  of 
current  sources,  each  of  which  is  determined  by  the  number  of  electrons  and  the  electric 
field  at  that  point. 


Figure  3.2  Generalized  equivalent  circuit  of  the  simulated  device  in  the  TOF  program, 
corresponding  tc  the  cross  section  of  the  MOSFET  in  Fig.  2.1  and  the 
MESFET  in  Fig.  2.12. 


3.1.2  Basic  Equations 

Two  basic  equations  describe  carrier  transport  (continuity  equation)  and  space-charge 
balance  (Poisson's  equation)  in  an  n-channel  FET  device.  Without  the  electron-hole 
generation  and  recombination,  the  two  equations  are: 


3nv 

IT 


(3.1a) 

(3.1b) 


with 


Jn  =  -  qov  +  qD(V\|/^j^)  Vnv 


(3.2) 
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where  ny  is  the  electron  density  per  unit  volume,  Vd,  is  the  electrostatic  potential  in  the 
channel,  referenced  to  the  Fermi  level  in  the  substrate  as  indicated  in  Fig.  3.1(c),  p  is  the 
net  charge  density  in  the  channel,  vfVyd,)  and  D(V\j/ch)  are  the  electron  velocity  and 
diffusivity,  respectively,  both  arc  function  of  electric  field.  Using  assumption  (1)  and  (3), 
and  noting  that  the  channel  is  in  the  y  direction,  the  above  equations  become: 

3n„  1  ^Nv 


dliy  _  2  ^Ny 

dt  q  3y 
p 


(3.3a) 


(3.3b) 


. _ ../^chx  .  dn, 

Jnv  *  q"v  v(-3 — )  +  qD(-3 — )  -3— 

Ny  ay  dy  dy  (3  4) 

Since  we  work  with  the  electron  density  per  unit  width  ni  on  each  grid  point,  the  final  form 
of  the  continuity  equation  after  discretization  is: 

n,(i,  t+At)  =  n  (i,t)  +  ( fy(i- 1 ,  t)  -  fy(i,t) ) 

W(i)  (35) 

where  i  is  the  index  for  the  discretized  position  along  the  channel,  and  At  is  the  self- 
adjustable  time  interval.  For  convenience,  the  electron  flux  in  the  channel  direction  fy(#/s) 
is  defined  to  be  in  the  opposite  direction  of  the  current  densities  jNy 

m 


fy(i,t)=-W(i)5- 


=  W(i)k(i.t)-5^ 


(n,(i+l,t)  -  ^.^(i,t))  0(14) 
Ay  Ay 


Basically,  the  program  solves  equation  (3.3)  -  (3.6)  at  various  instants  of  time  during  a 
cycle  of  oscillation.  Given  the  instantaneous  distribution  of  the  electron  density  and 
terminal  boundary  conditions,  the  program  determines  the  channel  potential  by  solving 
(3.3b).  Knowing  the  channel  potential  and  the  electron  concentration,  the  instantaneous 
electron  particle  current  (including  the  dependence  of  carrier  mobility  on  electric  field)  can 
be  calculated  using  (3.4)  from  which  the  time  derivative  of  the  carrier  concentration  can  be 
calculated.  From  the  time  derivative  of  the  carrier  concentration,  the  program  computes  the 
carrier  distributions  an  instant  in  time  later  (3.5),  and  repeats  the  cycle. 
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3.1.3  Algorithm 

The  algorithm  used  to  solve  the  above  equations  is  a  simple  forward  difference  scheme, 
i.e.,  the  electron  densities  at  time  t+At  are  calculated  in  terms  of  only  quantities  at  time  t  as 
indicated  by  Equation  (3.5).  The  coupled  equations  (3.5)  and  (3.6),  controlling  the  system 
variables  electron  concentration  and  the  flux  in  the  channel  directions,  are  updated  based  on 
their  values  of  previous  time  and  treated  separately.  This  requires  the  two  equations  be 
coupled  closely  by  small  updates  of  each  system  variable  [18]. 

Given  the  bias  Vsg  and  Vdg»  the  electron  densities  at  source  and  drain  can  be 
calculated.  Electron  densities  (#/cm)  at  other  nodes  ni(i)  can  be  assumed  to  be  the  linear 
taper  from  the  source  to  drain  at  time  t  =  0.  The  gate  voltage  VG(i)  is  ideal,  and  always 
varies  linearly  from  gatel  to  gate2  electrodes.  The  potential  in  the  channel  is  a  function  of 
the  electron  density  and  the  gate  voltage  as  determined  by  ID  Poisson’s  equation  (3.3b). 
The  electric  field  is  the  gradient  of  the  potential.  From  the  empirical  velocity-field  curve  and 
diffusivity-field  curve,  the  flux  in  the  channel  fy  can  be  obtained  through  Eqns  (3.4)  and 
(3.6).  At  this  point,  a  time  step  At  is  taken,  the  electron  densities  are  updated  according  to 
Eqn  (3.5).  If  the  electron  concentration  changes  more  than  0.5%  for  any  node  in  the 
channel,  the  time  step  is  reduced  and  the  electron  concentration  is  recalculated.  The  new 
flux  can  be  calculated  once  the  electron  density  is  accepted.  If  the  newly  calculated  flux 
diffen  from  the  old  flux  by  greater  than  1%,  the  entire  electron  concentration,  potential  and 
the  flux  calculation  are  redone  with  a  smaller  time  step.  This  keeps  the  currents 
(proportional  to  flux)  from  changing  the  electron  concentration  too  quickly  and  helps  to 
control  the  simulation.  One  cycle  of  the  time  iteration  is  done  when  the  flux  is  accepted. 

There  are  three  types  of  boundary  conditions  for  this  ID  problem.  One  can  either  (1)  fix 
the  electron  density  at  the  source  and  drain,  (corresponding  to  giving  Vsg  Vdg  to  find 
drain  current  Id),  or  (2)  fix  the  source  electron  density  and  drain  current  (corresponding  to 
giving  Vsg  and  Id  to  find  Vd),  or  (3)  fix  the  drain  electron  density  and  drain  current 
(corresponding  to  giving  Id  and  Vd  to  find  Vsg)^ 

The  flow  chart  of  the  TOF  program  is  shown  in  Fig.  3.3.  The  key  steps  are  discussed 
in  depth  in  the  following  sections. 


infliBtJciTsji 
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3.1.4  Calculation  of  Electron  Flux  Along  the  Channel 

The  flux  at  node  i  in  the  channel  direction  is  fy(i),  which  consists  of  a  diffusion  term 
and  a  drift  terra. 


where 


Ui)  =  fy‘‘'^(i)  +  f,‘^(i) 


Oi)  =  D(i)W(i) 


n,(i)  -  n,(i+l) 


for  V  >  0,  or 


.drift... 

fy  (l)=-^W(l)v(l) 


-frift...  ... 

fy  (0  =  -  --  W(l+ 1 )  V(l) 


(3.9a) 


(3.9b) 

for  V  <  0,  where  v  is  the  electron  velocity  along  the  channel,  which  is  in  the  opposite 
direction  of  the  channel  electric  field  ^ 

The  diffusion  coefficient  D(i)  and  individual  electron  velocity  v(i)  are,  in  general, 
nonlinear  function  of  both  electric  field  components  in  and  perpendicular  to  the  channel 
direction.  In  the  TOF  program,  only  the  tangential  electric  field  ^y  dependence  is 
considered.  The  drift  velocity  is: 


v(i)  =  •  Ho  5y(i) 


(3- 10a) 


(3.10b) 


v(i)  = 


Ho '  ^y(‘) 


L  ^  /  j  (3.10c) 

for  the  constant  mobility  case,  the  velocity-field  curve  of  Si,  and  the  velocity-field  curve  of 
GaAs.  Here,  %c.  is  the  saturation  field  of  the  velocity,  and  vjat  is  a  field  dependent 
coefficient  expressed  in  terms  of  fitting  parameters  Ay,  By,  Cy,  Dy  and  Fy: 
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Vs„  =  A,e 


1  + 


^o'^y 


The  diffusion  coefficient  is: 


(3.11) 


r  . 


(3.12) 


D(i)  =  Dg  +  ADe 

The  electric  field  in  the  channel  is 

e  /  X  ^(y)  _  V(i+1)-V(i) 

8y  -  Ay  (3.13) 

where  \}r(y)  consists  of  the  channel  potential  \|rch  and  a  non-local  potential  correction  Ax/ch 
due  to  the  influence  of  line  charge  at  neighboring  locations  in  the  channel.  For  the 
MOSFET  structure,  Vch  is  the  surface  potential,  while  it  is  the  potential  of  the  buried 
channel  for  the  MESFET  structure. 


3.1.5  Electrostatics 

When  the  gradual  channel  approximation  (GCA)  is  employed,  one  can  find  the 
electrostatic  potential  by  either  solving  the  ID  Poisson’s  equation  as  done  in  [19]  or 
working  with  Green's  function  [20].  With  GCA  removed,  the  Green’s  function  also  gives 
the  non-local  potential  correction.  In  Appendix  A  -  D,  the  later  method  is  used  to  derive  the 
channel  potential  as  well  as  the  non-local  potential  correction  for  the  RG-MOSFET  and 
RG-MESFET. 


3. 1.5.1  RG-MOSFET 

The  RG-MOSFET  is  a  surface  channel  device.  The  channel  potential  and  its  non-local 
correction  can  be  found  by  Green’s  function  (Appendix  B,  equations  (B-6)  and  (B-9)): 


Vd.(y)  =  V^y)  .3^+  V„  -  /vJ*2V,(Vo(y)-3^) 
'-’OX  V  '-'OX 


(3.14) 


=>>  27t(e,+e, 


,  /  .x2  2e, 

In(y-y)- 


msl 


Es+Eox 


•) . In[(y-y’)‘+4m24] 


( n,(y’)  -  n,(y) )  dy’ 


(3.15) 


V*’)  =  Vo(y).V„ 
Va(y)-Vo.+v<„i 

cX 


(3.16a) 

(3.16b) 


'-o*  (3.16c) 

AVch  is  due  to  the  non-uniform  distribution  of  surface  charge  in  a  real  MOSFET 
device.  After  discretization:  y -jAy,  y=iAy,  npnsAy,  and  summing  k  line  charges  on  both 
sides  of  the  point  of  interest, 


*  j«i-k 

ln((j-i)Ay)^-^^^(^ 


I  ln((j-i)Ay)‘--::^  >  (^^)"'-'ln(0--i)Ay)^+(2mt„)^) 

V  )  (3.17) 

In  the  simulation,  k  and  M  are  picked  such  that  any  further  increase  in  k  or  M  does  not 
increase  A\{rch(i)  substantially.  In  the  simulation,  k  is  around  7  and  M  is  around  6. 


3. 1.5.2  RG-MESFET 

The  RG-MESFET  is  a  buried  channel  device.  The  channel  potential  and  its  non-local 
correction  can  also  be  found  by  Green's  function  (Appendix  C,  equation  (C-10)  and  (C- 
11))  assuming  the  substrate  doping  Na=0: 

*  (3.18) 

Av^v(x,y)  =  qNp Jdy’  |  G(x,y,x’.y’)  dx’ 


where 


-4W  # 

J  G(x,y,x’,y')  dx'  =  [  (x-Xj(y))  In  ((y-y’)^+  (x-Xj(y))^ 


-  (x-Xj(y'))  In  ((y-y')^+  (x-Xj(y’))^ 

-  (x+x^(y’))  In  ((y-y’)^-t-  (x+Xj(y'))^ 
+  (x+x^(y))  In  ((y-y’)^+  (x+Xj(y))^ 


(3-19) 


45 


,  x-x  .(y)  ,  x-x  .(y’) 

+  2  ly-y-l  ( tan"'— , 
iy-yl  ly-yl 

_iX+x  (y’)  .i^c+x  (y) 

-  tan  — r-2-r-+  tan  i— )  ] 
ly-yl  ly-yl 


(3.20) 


Tqji  is  the  ion  implantation  depth  of  the  n-channel,  Nd  the  chainel  doping  density,  and  ns 
the  electron  density  per  unit  area.  If  a  non-unifonn  charge  distribution  (caused  only  by  the 
variation  of  the  depletion  region  width,  since  we  assume  the  channel  is  neutral)  is 
considered,  however,  the  resulting  channel  potential  varies  with  the  depth  of  the  channel, 
which  contradicts  the  original  assumption  that  channel  potential  is  only  a  function  of  y.  Fig. 
3.4  plots  the  channel  potential  correction  calculated  using  Eqn  (3.19)  and  (3.20)  for  two 
different  channel  depths:  the  depletion  edge  xj  and  the  centroid  of  the  channel  xm. 
(L=10jj.m,  Vgg=SV.  The  electron  density  in  the  channel  has  an  average  value  of  3.86x10^^ 
cm'2,  oscillation  amplitude  of  1.765x10^^  cm*2.  and  spatial  period  of  1.44}im.  The 
uncorrected  potential  varies  from  0.8V  to  8.8V).  Fortunately,  the  difference  is  small. 


Qiannel  Positon  (jj.m) 


Figiur  3.4  Channel  potential  correction  A\}rch  for  MES  structure  calculated  using  Eqn 
(3.19)  and  (  J.20)  for  two  different  channel  depths:  x=xd(edge  of  the 
surface  depletion)  and  x=Xtn  (centroid  of  the  channel). 
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For  further  simplification,  assume  that  the  channel  charge  is  a  5-function  at  the  centroid 

2K 


(3.21) 


The  corresponding  charge  distribution  is  shown  in  Fig.  3.5(a).  The  non-local  potential 

correction  can  be  calculated  with  the  aid  of  Fig.  3.5(b): 
j*i+k 

*  j»i-k 


j«i+k 

((j-i 


2  ni(j)*n,(i)  2 


2  n,(j)+n,(i)  2 


(3.22) 


metal 


Nd 

channel 


Na 
substrate 


(b) 


Figure  3.5  (a)  Charge  distribution  in  a  MESFET  assuming  the  channel  charge  is  a  5- 

function  at  the  centroid,  (b)  Calculation  of  non-local  potential  correction. 
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Fig.  3.6  shows  the  non-local  potential  correction  for  the  channel  centroid,  calculated 
from  Eqn(3.22).  Comparing  Fig  3.4  and  Fig.  3.6,  one  can  sec  that  the  two  methods  give 
approximately  the  same  results.  Mathematically,  the  non-local  potential  correction  using  the 
5  -  charge  distribution  at  the  centroid  is  simpler. 


Figure  3.6  Channel  potential  correction  Aych  for  the  MES  structure  calculated  at  x=xm 
(centroid  of  the  channel),  using  Green's  function  approach  (3.19)  &  (3.20), 
and  5-shcet  charge  approximation  (3.22). 

3.2  Simulation  Results  for  MOSFET  CDO 

The  MOSFET  CDO  simulation  has  been  performed  for  a  range  of  conditions  [21]. 
Type  (1)  boundary  condition  is  used,  i.e.,  the  electron  densities  at  both  source  and  drain 
arc  fixed.  The  nonlinear  velocity  field  curve  expressed  in  (3.10c)  is  used  with  the  following 
parameters:  no=5200  cm^/V-s,  Av=9. 55x10“^,  Bv=1.316xl0'7,  Cv=5.8xl0'^, 
Dv=6.788x105,  Fv=0.316.  The  constant  diffusivity  determined  by  the  Einstein  ralationship 
is  used. 

Fig.  3.7  shows  the  frequency  and  normalized  peak-to-peak  electron  density  as  a 
function  of  (nsd)'^  ns  is  the  average  number  of  electrons  per  unit  gate  area  in  the  channel. 
For  the  MOSFET  structure,  the  conduction  channel  is  at  a  constant  distance  d=tox  from  the 
resistive  gate.  This  relationship  implies  that  if  nj  is  halved  and  the  d  is  doubled,  the 
frequency  will  not  change. 
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It  is  expected  that  the  lower  frequency  limit  is  determined  by  the  maximum  charge 
capacity  of  the  potential  well  in  the  devices.  For  an  ideal  MOSFET,  the  inversion  electro*' 
density  can  be  very  high.  However,  to  implement  a  MOSFET  structure  in  reality  using 
GaAs  semiconductor,  an  AlGaAs/GaAs  heterojuction  has  to  b«  used.  The  energy  barrier 
provided  by  band  discontinuities  is  much  smaller  than  that  of  a  ideal  insulator.  The  leakage 
current  across  the  device  (along  the  x  direction)  due  to  thermionic  emission,  which  is 
inversely  proportional  to  the  exponential  of  the  barrier  height,  will  be  much  greater.  As  a 
result,  a  high  inversion  electron  concentration  in  the  channel  is  more  difficult  to  obtain. 

As  indicated  in  reference  [22],  the  maximum  2DEG  density  is  determined  by  the  doping 
density  of  AlGaAs  and  GaAs  and  by  the  conduction  band  discontinuity  at  the 
AlGaAs/GaAs  interface.  It  is  limited  in  the  ncighboriiood  of  10^2  cm*2  for  AlGaAs  doping 
ranging  from  5x1017  to  IQl^cm'^.  The  thickness  of  the  AlGaAs  is  chosen  so  that  the 
surface  depletion  region  overlaps  the  interface  depletion  region  to  allow  the  gate  modulation 
of  the  channel  so  that  there  are  no  free  carriers  in  the  barrier  layer  and  the  channel  is 
modulated  by  the  gate. 

The  minimum  frequency  is  then  determined  by  the  maximum  surface  electron  density  in 
the  channel  (-10l2cm’2)  and  the  thickness  of  AlGaAs  (~800A),  which  gives  l/(nd)~12.5A, 
corresponding  to  a  frequency  of  about  50GHz. 

The  maximum  frequency  can  be  obtained  by  examining  the  rate  of  change  of  the  charge 
packet  width  [23]: 

-  3nj  ,  [I  n  . 

<3.23) 

where  D  is  the  diffusion  constant,  ^  tne  effective  differential  mobility,  and  C  the 
capacitance,  y  is  along  the  channel  direction.  The  first  term  in  (3.23)  represents  the  charge 
packet  expansion  due  to  the  electron  drifting  under  the  influence  of  the  self-induced  electric 
field.  The  second  term  represents  the  charge  packet  expansion  due  to  the  electron  diffusion, 
which  is  always  positive.  When  the  differential  mobility  is  negative  (^<0),  Wp  could  be 
negative,  indicating  the  charge  packet  contraction  instead  of  expansion.  i.e.. 


*  *  (3.23) 

For  typical  values  of  D=135  cm^/s  (D=kBT/q  po.  ixo=5200  cm^/V-s),  ~  -850  cm^A^-s  at 
a  field  of  approximately  1  V/pm,  l/(nsd)max  ~  88  A.  From  Fig  3.7,  however,  the  oscillation 
amplitude  decreases  dramatically  as  the  oscillation  frequency  increases.  Therefore,  the 
practical  maximum  frequency  limit  is  determined  by  the  oscillation  amplitude. 


0  »  iO  JS  30  55  30 

(fi'ir*  lA*) 


( ’iiiviuisticn  irstilu  lot  ilif*  MOM)  I  (  I  m  )  i  i  lu  v 

pc.ik  t(>  jir  ik  ifl'plifii.Ir  a  (iilH  Uori  o(  '!ir  pn  jiin'Vf  (n  .! 

il  liif Niiip  ih.ii  ih^  (■  niil  n  1  0  ^ 

IS  '• !)  1 1  V  c'ii  ifri  jl'K’IU  V  !  .* !  t 


Mil  in  t-'r  MiiU  '•  >r  t'iT  \1l  '« )  I  I  (  I  X  ) 

'.!.M  jrsii.'n  iiiutt  il,r  Nfl  t  f  (i.iiviiro!  ud  sSm  sn  m  i  i,;  i  ■; 

i  ,,,  V  I!  !i  ii)'(  ili'pinj;  liriuitir  s  .ur  ms jsc.!  ) ''i- 

li  if  VC,'. «  i!\  licl.l  I.  ufvr  .(tr  :!if  v.jitx”  ,u  UNf'ii  m  stif  lit  I  M<  .1  )  I  .  .i  :,; 

:  IS*  i  '.I'i  ..  .litC!  iO  ’>!<'  il.Ua.  O't*  i .11  H  .()  (  M' .  j '  if  Ht  V  l  <  .t  | '1 'l  ■  1  *,  H  S  :.1  [f  f  •.  I 


50 


function  of  (nsd)  ’  up  to  26A  and  then  saturates  thereafter.  Here  d=x<i  is  the  separation 
distance  measured  from  the  edge  of  the  depletion  region  to  the  resistive  gate/GaAs 
interface. 
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(a)  l'tr<)ufncy  as  fiim  tion  of  the  parameter  (nd)  •  (A), 
th)  ( )>.<  illation  afnpiiiude  as  function  of  n,  (cm  •') 


Assuming  that  the  charge  goes  into  the  channel  as  a  rectangle  with  height  equal  to  the 
channel  doping,  the  depletion  width  moves  towards  the  gate-GaAs  interface  as  the  channel 
gets  fuller  and  away  from  the  interface  as  n  gets  smaller.  Therefore,  xj  oscillates  as  electron 
density  oscillation  occurs,  and  d  is  no  longer  a  constant  as  that  in  MOSFET  device.  The 
computer  progr.am  utilized  the  assumption  that  even  though  the  electrons  can  redistribute 
inside  channel  (along  the  x  direction),  the  flux  or  cunent  is  still  one  dimensional  (along  the 
y  direction). 

It  seems  that  the  relation  is  independent  of  the  doping  density  and  it  is  obvious  that  the 
frequency  can  be  modulated  in  a  wide  range  by  varying  the  average  sheet  electron  density 
n^,  w  hich  is  controlled  by  the  source  to  gate  voltage  VsG-  A  high  frequency  of  about  1 10 
GIIz  is  achievable  for  the  MESFET  CDO. 

The  oscillation  amplitude  l/2Anpp  is  approximately  a  linear  function  of  sheet  electron 
density  n^.  As  the  frequency  increases,  the  a.c.  output  power  decreases. 

Tlie  minimum  frequency  fmin  is  dctcmiined  by  the  maximum  njd  product.  Since 


^  D  (3.24) 

the  maximum  value  of  n^d  (obt.aincd  by  differentiating  (3.24)  with  respect  to  n?)  occurs 
w  hen  the  channel  is  "half  full"  (neglect  the  zero  bias  depletion  width  Xtio): 


(3.25) 


(3.26) 


f  - 

mm  IM  p 

Ti  (3.27) 

where  A  is  the  slope  of  the  f  vs.  (n^d)'*  plot  for  l/(nx,i)  <  25A.  Therefore,  fmm  i'^ 
dctrnnincd  by  the  l.ihncauon  parameters,  NoandTepi. 

Ilu*  maximum  lieqi.ency  is  limited  to  1  KXjllz.  Tlie  reason  is  not  yet  known. 
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3.4  Summary 

The  one  dimensional  Time-Of-Flight  program  is  used  to  simulate  the  resistive  gate  FET 
devices.  By  combining  the  transferred  electron  effect  (nonlinear  velocity  -  field  curve)  in 
GaAs  with  the  two  dimensional  geometry  of  a  FET  device,  a  series  of  contiguous  electron 
packets  can  be  generated  in  a  resistive  gate  FET.  High  frequency  oscillation  in  the  GHz 
range  can  be  obtained.  The  frequency  and  the  oscillation  amplitude  depend  on  the  channel 
electron  density  and  the  distance  between  the  gate  and  the  channel.  The  electron  density  can 
be  changed  by  changing  the  source-to-gatel  voltage. 
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CHAPTER  4 

EXPERIMENTS  ON  THE  RESISTIVE  GATE  MESFET  OSCILLATOR 


The  resistive-gate  MESFET  structure  was  first  fabricated  in  order  to  verify  the 
contiguous  domain  oscillation  concept  since  it  is  the  most  widely  used  technology  for  GaAs 
integrated  circuits  in  industry.  The  devices  were  fabricated  at  ITT  (International  Telephone 
and  Telegraph)  Gallium  Arsenide  Technology  Center  based  on  the  design  and  process 
dt  'eloped  here  at  Purdue.  Three  wafers  with  exactly  the  same  structure  but  different 
channel  doping  and  channel  depth  were  fabricated  (CDO-1-3:  channel  doping  =  2.4x10*^ 
cm‘3,  depth  =  20(X)A;  CDO-1-4:  channel  doping  =  tixlO'^cm*^,  depth  =  ISOOA;  CDO-1-5: 
channel  doping  =  1.2x10^^  cm'3,  depth  =*  ISOOA.  Source  /  drain,  doping  =  SxlO*"^  cm'3, 
depth  =  2(X)0A  for  all  the  samples).  Only  wafer  CDO-1-5,  which  exhibits  microwave 
oscillation,  is  considered  here.  The  other  two  wafers  did  not  oscillate.  The  information  on 
the  other  two  wafers  of  device  can  be  found  in  Ref.  [19]. 

This  chapter  summarizes  the  experiments  on  the  Resistive-Gate  MESFET  devices.  The 
fabrication  sequence  of  the  RG-MESFET  (CDO-1-5)  is  covered,  followed  by  dc  and  the 
microwave  testing  results  and  discussion.  Further  modification  of  the  ITT  devices  is  also 
presented. 


4.1  Fabrication  F^rocedure  of  the  RG-MESFET  CDO 

The  device  structure  of  a  fabricated  resistive-gate  MESFET  oscillator  is  shown  in  Fig. 
4.1.  The  basic  device  dimensions  are  the  channel  depth  Tepi,  gate  length  Lq,  gate  width 
Wg,  channel  length  L,  channel  width  W,  and  the  overlap  of  the  resistive  gate  over  the 
source  or  drain  region  X.  The  origin  of  X  is  at  the  source  (drain)  -  channel  junction.  When 
the  ohmic  contact  is  on  top  of  the  source  (drain)  region,  X  is  positive.  If  the  ohmic  contact 
is  on  top  of  the  channel,  X  is  negative. 

The  complete  fabrication  of  the  RG-MESFET  device  requires  eight  mask  levels, 
summarized  in  Table  4.1.  The  process  starts  with  a  two-inch  semi-insulating  GaAs  wafer 
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with  (100)  orientation.  A  350A  PECVD  SiOH  layer  is  first  deposited  to  minimize  the 
channeling  effect  during  the  ion  implantation.  The  first  lithographic  step  is  the  registration 
mask  (REG  mask),  which  defines  alignment  marks  for  the  subsequent  masks.  Since  the 
first  two  steps  of  the  fabrication  are  ion  implantation  for  the  channel  and  the  source/drain, 
which  leave  no  visible  marks  on  the  GaAs  wafer,  it  is  necessary  to  have  this  registration 
mask  to  etch  the  GaAs  before  the  ion  implantation.  CF4  Plasma  (200mTorr,  150W,  6mins) 
was  used  to  etch  the  SiOH  cap  and  expose  the  GaAs  that  was  to  be  subsequently  wet 
etched  to  form  the  registration  marks. 


Figure  4. 1  Schematic  cross  section  of  the  RG-MESFET  Device. 

I 


Table  4. 1  RG-MESFET  device  p.'Dcess  summary. 


Operation 

Mask  level  (field) 

1.  350  A  PECVD  SiON  ion  implant  cap  dcp>osition 

2.  Registration  mark  etch 

REG  (dark) 

3.  Selective  channel  implant 

Ln(dark) 

4.  Selective  source  /  drain  implant 

NPI  (dark) 

5.  Anneal  cap  deposition,  anneal,  cap  removal 

6.  1000  A  PECVD  SiN  deposition,  patterning  and  etch 

SIN  (clear) 

7.  Au/Ge/Ni  ohmic  contact  evaporation,  liftoff  and  alloy 

OHM  (clear) 

8.  Resistive  gate  spuner  deposition  and  liftoff 

GAT  (clear) 

9.  Ti/Pd/Au  interconnect  evaporation  and  liftoff 

MTl  (clear) 

10.  2050A  PECVD  SiN  deposition,  pattemine  and  etch 

NVA  (dark) 
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The  channel  was  patterned  by  the  second  mask  (LFI  mask)  for  channel  ion 
implantation.  Doubly  ionized  silicon  Si++  was  first  implanted  into  the  GaAs  with 
acceleration  energy  of  155keV  and  dose  of  3.1  IxlO'^cm"^.  The  second  implantation  with  a 
lower  acceleration  energy  (65keV)  and  lower  dose  (5.2x10' ^cm'^)  was  then  done  to  retain 
the  carrier  concentration  toward  the  surface.  The  third  mask  (NPI  mask)  defines  the 
source/drain  implantation.  Si"'"'"  was  implanted  into  the  GaAs  with  an  energy  of  200keV 
and  dose  of  2xl0'3cm‘2.  Activation  is  subsequently  done  in  order  for  Si  at  interstitial  sites 
to  move  to  substitutional  sites  and  for  the  wafer  to  recover  from  the  implantation  damage. 
To  achieve  a  uniform  and  controllable  activation,  another  cap  layer  of  2000A  SiOH  is 
deposited  before  the  process  and  removed  with  5%  HF  for  15  mins  after  tlie  process.  The 
top  view  of  mask  levels  2  and  3  is  shown  in  Fig.  4.2,  along  with  the  cross  section  of  the 
device  with  the  channel  and  source/drain  ion  implantation  step  completed. 
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S.  I.  GaAs  Substrate 


Figure  4.2  (a)  Top  view  of  mask  1  (LFI)  and  mask  2  (NPI).  (b)  Cross  section  of  the 

^vicc  with  channel  and  source/drain  ion  implantation  completed. 


The  fourth  mask  (SIN  mask)  defines  a  ring  of  silicon  nitride  film  around  the  channel.  It 
reduces  the  field  enhancement  at  the  gate  edge  to  prevent  the  reduction  of  the  planar 
breakdown  voltage  of  the  gate-channel  Schottky  diode  [191.  The  silicon  nitride  is  LPCVD 
deposited  to  a  thickness  of  1 1(X)A.  The  unwanted  part  of  the  silicon  nitride  is  etched  away 
using  CF4  plasma  (4(X)mTorr.  6mins).  The  device  with  silicon  nitride  ring  formation 
completed  is  shown  in  Fig.  4.3. 


Nitride 


S.  I.  GaAs  Substrate 


Figure  4.3  (a)  Top  view  of  related  mask  levels,  (b)  Cross  section  of  the  device  with 

silicon  nitride  ring  formation  completed 


Mask  5  (OHM  mask)  defines  the  ohmic  contact  in  the  source  and  drain.  The  ohmic 
contact  is  formed  by  a  Ni/Ge/Au  (150/3()0/2()(X}A)  evaporation  and  heat  pulsed  at  390^0 
for  20  sec.  The  device  with  ohmic  contact  formation  completed  is  shown  in  Fig.  4.4. 

Mask  6  (GAT  mask)  defines  the  resistive  gate  pattern.  The  gate,  made  of  WSiN 
resistive  film,  is  deposited  by  rf  reactive  sputtering  from  pure  W  and  Si  targets  in  an  argon 
and  nitrogen  (65%  partial  pressure)  atmosphere  (power  =«  1.5kW,  pressure  =  10  mTorr, 
time  =  65  mins).  The  resulting  film  has  a  sheet  resistivity  of  30kQ/squaie  and  a  thickness 
of  1  lOOA.  The  Schottky  contacts  were  found  to  have  ideality  factor  of  1.65  and  a  barrier 
height  of  0.6V.  A  Lift-off  process  was  used  to  pattern  the  resistive  gate.  The  device  with 
resistive  gate  lift-off  step  completed  is  shown  in  Fig.  4.5. 
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Figure  4.4  (a)  Top  view  of  related  mask  levels,  (b)  Cross  section  of  the  device  with 

ohmic  contact  formation  completed. 


Figure  4.5  (a)  Top  view  of  related  mask  levels,  (b)  Cross  section  of  the  device  with 

resistive  gate  lift-off  step  completed. 
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Mask  7  (  MTl  mask)  defines  the  metal  connection  to  the  devices.  Ti/Pd/Au 
(500/500/5000A)  is  evaporated  and  patterned  by  lift-off.  The  device  with  Ti/Pd/Au  step 
completed  is  shown  in  Fig,  4.6. 
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Figure  4.6  (a)  Top  view  of  related  mask,  (b)  Cross  section  of  the  device  with  Au/Ge/Ni 

metalizadon  step  completed.  I 


The  last  layer  deposited  on  the  wafe”  is  a  cap  layer  oflLPCVD  SiN  to  protect  the 
devices  from  scratch  and  contamination.  Windows  are  pattemlwith  CF4  (400mTorr,  75W, 
llmins). 

There  arc  eighteen  devices  on  each  die  with  different  channel  length  and  width,  where 
L=10,  20,  50  fxm  and  W=10,  20,  50,  200,  and  400  jim.  The  offset  X=5  fim  on  both 
source  and  the  drain  sides.  The  layout  is  shown  in  Fig.  4.7. 
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Figure  4.7  The  layout  of  RG-MESFET  devices.  The  channel  length  and  width  LxW 
from  left  to  right  for  the  devices  on  the  top  row  are:  10x10,  20x10,  50x10, 
10x20,  20x20, 50x50.  LxW  for  the  devices  on  the  second  row  are:  20x400, 
10x200,  20x200.  LxW  for  the  device  on  the  third  row  are:  50x400, 
50x200,  20x200.  The  device  on  the  forth  row  are  the  same  as  that  on  the 
first  row.  Besides  eighteen  RG-MESFET  devices,  there  are  resistance  test 
patterns  for  gate,  channel  and  source/drain. 


4.2  rV  Measurement  of  the  RG-MESFET  Devices 

The  DC  electrical  testing  was  conducted  for  the  oscillator  devices  before  any  microwave 
testing.  A  Hewlet-Packard  4145B  semiconductor  parameter  analyzer  is  used  to  measure  the 
gate-to-source  and  gate-to-drain  Schottky  diode  characteristics  and  the  IV  characteristics  of 
the  RG-MESFETs.  The  drain  breakdown  voltage  is  found  to  be  around  4.5V.  Typical 
drain  current  Id  and  gate2  current  Ig2  vs.  Vgq  with  Vse  as  parameter  are  shown  in  Figure 
4.8.  Two  features  can  be  noticed  on  these  curves:  (1)  At  small  Vgg.  application  of  Vsg 
delays  the  turn  on  of  the  drain  current.  (2)  A  strong  negative  transconductance  (i.e.,  drain 
current  decreases  as  gate  field  increases)  is  observed  at  higher  Vgg*  This  behavior  can  be 


60 


understood  with  the  help  of  the  DC  equivalent  circuit  [25]  as  shown  in  Fig.  4.9  and  the 
schematic  potential  diagrams  for  various  bias  conditions  in  Fig.  4.10.  In  the  equivalent 
circuit,  Rx  denotes  the  gats  resistance  between  terminal  G1  and  point  x,  directly  over  the 
edge  of  the  source  N+  implant  region. 


Iq  (mA)  l02  (mA) 


Figure  4.8  The  drain  cuirent  Id  and  gate2  current  Ig2  vs.  Vcq  with  VsG  ^  parameter. 

Vdg  is  held  at  0  and  has  little  effect  on  the  drain  current  Id-  Vsg  increases 
from  0  to  6V  with  steps  of  IV. 
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■  Channel 


lD  =  qWnvV(^) 


Figure  4.9  DC  equivalent  circuit  for  the  RG-MESFET.  Rx  represents  the  gate 
resistance  between  terminal  G1  and  the  edge  of  the  source  [25]. 


Refering  to  Fig.  4.10(a)-(c),  where  non  zero  Vsg  is  fixed  at  a  value  less  than  the 
channel  pinch-off  voltage  Vp(ch),  the  drain  current  Id  is  initially  less  than  zero  (Fig 
4.10(a)),  since  Vdg  =  0-  Here 

is  defined  as  the  channel  pinch-off  voltage,  where 

kj,T  Nc 

^  (4.2) 

is  the  built  -  in  potential  of  the  Schottky  diode,  <j>Bn  is  the  barrier  height,  and  Nc  is  the 
effective  density  of  states  in  the  conduction  band.  With  the  application  of  Vgg.  however, 
the  drain  current  increases.  When  Vsg  =  Vdg  +  Vgg.  the  drain  current  is  zero  (Fig. 

4.10(b)).  The  gate  voltage  required  for  this  turn  on  of  positive  drain  current  is  simply: 

VQQ(on)  =  VjQ  -  VpQ  VjQ  <  Vp(ch)  . . 


Vso<Vp(ch) 


For  any  Vgg  greater  than  the  above  VGG(ott).  the  drain  current  is  positive  as  shown  in  Fig. 
4.10(c). 

Now  refer  to  Fig.  4.10(d)-(g),  where  Vsg  is  greater  than  Vp(ch).  The  channel  is 
pinched  off  at  the  source  and  the  drain  current  is  negative  when  Vgg  =  0  (Fig-  4.1 1(d)). 
The  application  of  a  gate  voltage  Vgg  increases  the  drain  voltage  as  well  as  the  potential  at 
point  X  due  to  the  resis’^ance  Rx-  When 


V  = _ S _ 

CG  L0-X3 


(Vp(ch)  -  V„) 


(4.4) 
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Ili'sh  i'tjii,»ii(i!',s  (4  3)  (4  It)  inil,i  :*tr  (iiat  as  V<,(>  is  nu'tf  ,ni  il,  ihr  i.jfn  <  o  sn’t.ii'-'  !"r 

V\,^;  is  iPi  rrasrd. 

As  IS  f  jrihcr  UHtrascd,  fn  nuirasrs,  and  esfiitn.dlv  Ix  ii'i'ir  mis f,  i  r 

dH'  Si  lu'iAy  die  If  at  tlic  sourrr  end  is  torvsan.l  (nasrd  Sint  r  i.  lianncl  n  d  iSc  i:  >’;•  n'  ''»* 
s  ir 'V  rd  as  a  disip'.nncil  Stduntk y  duHlr,  the  disttd’ii'ed  i. uirrnts  i  an  f1i >vv  !  -  i  .;n  ih<'  ,  liii'ri  ; 
to  d'C  y,a'r  nndnt  tins  uti  umstances  and  spdlosiT  oi  i  urs  I  ho  nu  r'“aii(ti'  i;  it;*  i  lini-ni  1, ,  j 
duo  to  vpijimrr  tt suits  in  the  doi.  tease  in  di am  t  nnrnt  Ip  The  no  aai, no  tr  an  a  I'n-diji.  t  m. 
ot  i()  is  ol'-.orvrd  Iho  high  lie  Id  negative  molnliiy  may  also  eontnhiite  to  (fio  ncaoivo 
iransoondiK tain 0,  but  ibis  elfoi  t  is  small  sime  not  nuu  h  nei-ative  tia.nsi  on  ia.  lan,  f  o, 
observed  (MK  e  the  do  vice  is  itKHhhed  to  set  Xyd) 

In  rig.  4  11,  the  lour  terminal  cunents  l.y,  Ij),  loi,  U,;  air  sli.nan  as  *  n.  'ions  'l  'd,> 
souive  to  gatcl  voliare  Vytj  andiitam  to  gate2  voltage  V|h;  As  dioHn.  •.I'lilo*.  cr  o.  luts 
at  die  soiin  c  when  V'sii  <  d.'tV.  whiih  is  maiulcsied  by  a  laignr  I,  ilnn  Ip  In  tho  i.  a 
the  .vile  to  source  diode  is  effectively  fotwaird  hiascd.  With  the  imirise  ol  \  ,j,.  the 
ei!-itive  bias  voltage  across  the  soiiri,'e/<  hannel  dicule  V'xy  is  n'dut  ed  and  the  dnsle  o 
ini.froe  biased  eventually,  1  urther  increase  of  Vm;  causes  I^  and  l|,  to  dri.!'-t:,r 
-lenti.iHv,  As  V\s  approaches  the  hirakdown  voliage  Vnu  whu  h  is  about  fi  .'A  tot 
\l,  •  A  10'’  I  m  ’  10|,  (he  generation  of  elec  tion  hole  pairs  through  as  al  u'l  be  bp'  iPluvc  u 

Hl’llfS 

.At  die  drain  riicl  (l  ig.  4, 1 1  (h)).  Ip  t>eIow  4  3 V  is  weakly  d»'i'endent  on  Vi «  shoe*  m,* 
that  the  device  can  iW  approximated  as  a  current  source  in  tins  range  1  or  V’lx;  c'caier  iluin 
4  3V,  tlie  drain  cument  increases  rapidly  and  the  gaie2  cuoent  decieases.  ind.i  attne 
avalanc  he  hrrakdowii  of  the  gate  to  dram  junction  Hie  lower  V'tx;  value  where  btrak.ioccn 
(H  curs  near  ilie  drainUalMuit  4V)  than  that  of  Vscj  where  breakdown  cvviiis  near  tbe  souii  e 
(al'out  7V)  is  caused  iW  the  gate  ovei  lapping  of  drain  rriMon  (late  oveilappine  lu  tlie  di  un 
trgion  etfeciively  biasits  the  gate  to  dram  diode  to  a  higher  value  than  tlie  !|'pheil  Vi„ ; 

1  tom  the  above  discussion,  it  is  shown  that  the  eKiension  of  the  rr  ustive  potoon  of  the 
gate  into  the  source  region  has  several  disadvantages.  I'irst,  a  huger  Vc;  need.s  to  he 
applied  as  comp.nrrd  with  the  same  LxW  deviiT  with  Xs”Xi)-d)  m  order  to  achiece  the 
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r!(\  lilt'  (ickl,  'iincr  ihf  g-itc  is  lonj^rr  Second,  the  gate  voltage  tends  to  forward  bias 
itie  4atc  ii>  si'on  c  dusir  due  to  the  v<>lt,ige  dtop  on  Rx  ITiCtefote,  a  larger  Vsg  >*  required 
to  rroiiuT  the  duxle  is  reverse  biased  Spillover  of  clectnins  fremi  thannel  tr>  gate  van  ixcur 
when  (he  ^ate  t('  source  Scboltky  rSiixlc  is  forward  biased  Lastly,  the  effective  (internal) 
IS  voitj'lrd  to  Ihf  s<Hin  f  election  density  is  ilcicrrnincd  not  only  t’v  the  esrcrnal 
S'SC  but  also  by  V(ir;, 
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figure  LI  I  Llfcct  of  source  to  gate  1  voltage  V<;r,  and  dram  to  gatc2  voltage  VpG  on 
four  tcnninal  curTenis  for  .i  10x10  (LxW)  (T){^-l-5  device.  A  Vgo^HiV  is 
applied,  whii  h  corresponds  to  an  electric  field  of  S(XX)V/cm. 

(iO  Is.  I(il.  If).  loj  vs,  Vsrj. 

(Id  Is.  I(ii.  Id.  I(:2  vs.  V,k;. 

4  d  Microwave  rharactcil/aiion  of  the  RG-VlLSniT 
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I  'lc  w  hrin.itii  bias  circuit  setup  is  shown  in  big.  4.12.  The  saniiiic  is  mounted  on  i 
!  ( can  usuu:  clcctnc  conductive  epoxy  or  silver  paint  to  ensure  a  gcxxi  contact.  All  the 
connci.  tions  ue  made  by  aluniinutn  wires.  Idic  source,  jratcl,  and  gaic2  leads  arc  binutdu 
out  dic  I'ottorn  cl  r(.)5  Can  lor  bias,  while  the  dram  lead  diTctly  connects  to  a  striplinc  with 
t  hat  1.,  icMsnc  imped. nice  of  5(U2.  An  SMA  connector  couples  the  output  current  out  to  a 
bi  is  tee.  Hie  dc  output  is  connected  to  the  DC  biasing  circuit,  and  the  ac  output  lead  of  the 
bias  lee  i:;  ctsfuiected  to  the  s(vctrum  nnalyrer  by  coaxial  cable  for  frrcpicncy  measurement 
Ik-Iow  21  (dl/  (internal  mixer  is  used).  I'or  fiequencies  higher  than  21GH/.  an  external 
«  :uf guide  mixer  is  used  via  a  lunchcr  which  connects  the  signal  from  the  striplinc  to  tlic 
vs  ueiMiide 


I  igure  t  12.  Schematic  diagram  for  bias  cia'uit  setup. 

W'ficn  switch  S  is  at  "1",  the  signal  is  directiv  connected  to  the  internal 
mixer  of  the  spectrum  analyzer.  When  .S  is  at  '2'’,  a  luneber  and  an  extenuil 
waveguide  mixer  are  used  between  he  signal  and  spectrum  analyzer. 


To  adties  e  successful  operation,  projicr  bias  con  tion  is  required.  First  of  all,  the  g.ite 
voltage  V(;(;  has  to  bias  the  channel  beyond  the  threshold  for  the  electnm  transfer  effect. 
Volt.ige  Vse,  should  be  such  that  it  reverse  biases  the  source  and  gate  junction  to  prevent 
the  spilliiver  of  electrons  from  channel  to  gate.  Tlic  higher  the  reverse  hia.s  voltage  V'sd  is 
.u  toss  the  source  gate  1  terminals,  the  fewer  electrons  will  l>c  injected.  Die  dram  and  g.ite  2 
junction  also  needs  to  be  reverse  biased  in  order  to  isolate  the  channel  and  the  gate  as  wcil 
as  the  channel  and  the  dram. 


Fig.  4.13  is  the  pwwcr  spectrum  for  a  10x10  device.  The  fundamental  frequency  is 
found  to  be  12.933  GHz.  Oscilladon  signals  arc  also  found  at  25.899  GHz  and  38.857 
GHz.  which  arc  the  first  and  second  harmonics. 


Figure  4.13  Power  spectrum  for  a  10x10  itm  RG-MESFET  device  shown  in  Fig.  4.1. 
VsG"5.5V,  Vgg“I6V,  Vdg=^V.  The  overlap  Xs*XD=5iim. 


The  dependence  of  oscillation  frequency  and  power  on  bias  voltage  Vgg,  Vsg,  and 
Vdg  arc  .shown  in  Fig.  4.14.  It  is  obvious  that  both  Vgg  and  Vgo  have  strong  effects  on 
the  device  oscillation,  while  Vdg  has  little  effect. 
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Vgg  can  affect  the  operation  frequency  in  two  ways.  First,  any  chang.  in  Vgg  will 
change  the  electric  field  along  the  gate,  which  changes  the  drift  velocity  of  eU  /ns  in  the 
channel.  Second,  since  the  resistive  gate  overlaps  the  source,  changing  Vgg  'auses  the 
internal  gatc-to-source  potential  to  change  in  the  same  direction. 

The  oscillation  frequency  varies  with  Vsg.  from  a  minimum  of  12.7  GHz  at  5V  to  14.7 
GHz  at  6.1V.  The  frequency  tunability  is  about  2GHzA^. 

Vdg  has  very  little  effect  on  the  oscillation,  indicating  the  device  can  be  viewed  as  a 
constant  microwave  source,  consistent  with  the  current  measurement  results  in  Fig.  4.1 1. 
As  shown  in  Fig.  4.14,  the  output  power  of  a  10-)im-wide  device  ranges  between  0.1  and 
1  |4W  for  frequencies  around  13  GHz.  This  low  power  is  purely  due  to  an  exponential 
decay  of  the  ac  current  with  distance  under  the  large  gate2  contact. 

The  fundamental  frequency  is  in  the  range  of  10-20  GHz,  which  could  be  explained  by 
a  Gunn  diode  of  10  ^tm  channel  length.  In  order  to  identify  the  oscillation  mode,  the 
frequency  dependence  on  the  channel  length  is  investigated.  The  oscillation  frequency 
versus  mverse  channel  length  for  four  devices  is  shown  in  Fig.  4.15.  The  frequency 
increases  linearly  with  inverse  length,  having  a  slope  of  1.3x10"^  cm/s.  The  solid  line  is 
calculated  from  f  =»  vj  /  L,  representing  a  Gunn  diode.  The  scatter  of  the  data  is  caused  by 
many  factors  such  as  the  misalignment  of  the  gate,  different  injection  levels  at  the  source 
Vsg  2ind  different  applied  electric  field,  Wqq/L.  The  dependence  of  frequency  on  the 
channel  length  L  indicates  that  the  oscillation  is  due  to  a  single  domain  rather  than  multiple 
domain  formation. 


Figure  4.15  Frequency  vs.  (1/L)  plot. 


►  \ 
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4.4  Discussion 

4.4.1  Properties  of  tlie  RG-MESFET  Oscillator 

The  microwave  measurement  indicates  that  the  RG-MESFET  is  operated  in  a  single 
domain  mode.  Is  this  RG-MESFET  oscillator  simply  a  Gunn  oscillator?  Tlie  following 
observations  indicate  that  the  answer  is  no. 

Take  a  device  with  Xs  =  Xq-  0  and  LxW  =  20  x  20  )im.  Since  the  gate  overlap  of  the 
source  is  eliminated  ,  Vsg  ii^d  Vgg  arc  decoupled.  Variations  in  Vgg  will  not  affect  the 
electron  injection  level  at  the  source.  To  see  the  effect  of  Vqg  on  fhc  oscillation,  the  sourcc- 
to-gatel  voltage  Vsg  >s  fixed  at  0.5V.  Table  4.2  summarizes  the  observation  results. 

Notice  that  in  Table  4,2,  the  effective  drain-to-source  bias  Vds^Vgg+Vdg-Vsg  is 
maintained  at  13.5  V.  Oscillation  only  occurs  when  the  gate  field  Vgg  /  (L+Ks+Xd)  is  in 
the  regime  of  negative  differential  nxability  for  electrons  in  GaAs,  in  this  case,  greater  than 
10.5V  /  20  pm  =  5250  V/cm.  Therefore,  the  channel  field  is  not  detennined  by  the  dmin-to- 
source  voltage,  but  rather  by  the  linear  potential  drop  along  the  resistive  gate. 

For  the  RG-MESIEiT  oscillator,  the  frequency  can  be  tuned  by  Vsg.  although  the 
tunability  is  not  very  impressive.  The  frequency  can  be  changed  by  changing  the  electron 
injection  level.  The  fact  that  the  frequency  is  not  affected  by  changing  Vjx)  proves  that  the 
device  can  be  approximated  as  a  constant  current  source,  as  expected  for  the  CDO  device. 
The  only  function  of  the  drain  is  to  remove  any  electrons  which  drift  down  the  channel. 


Table  4.2  Effect  of  bias  condition  on  the  oscillation  behavior.  Device  parameters:  L 

X  W  =  20  X  20  pm,  Xs=Xd=0.  The  electron  injection  level  is  held  by 
applying  a  constant  source-to-gate  1  voltage  Vsg=^.5V.  Vgg  is  stepped, 
and  V[xj  is  adjusted  at  each  bias  to  keep  V^s  constant  at  13.5V, 


Vgg  00 

Vrxi  (V) 

Vds  (V) 

Oscillation? 

14 

0 

13.5 

Yes 

12 

2 

13.5 

Yes 

10.5 

3.5 

13.5 

Yes 

10 

4 

13.5 

No 

9 

5 

13.5 

No 

7 

7 

13.5 

No 

70 


4.4.2  Requirement  of  the  Resistive  Gate 

The  natural  questions  are  why  there  is  no  contiguous  domains  forming  with  the 
resistive  gate  on  top,  and  if  the  resistive  gate  is  really  screening  the  space  charge. 

In  the  computer  simulation,  the  linear  potential  on  the  gate  is  always  fixed.  Any 
influence  of  space  charge  effects  on  the  gate  potential  distribution  is  neglected.  If  this 
influence  is  calculated,  assuming  an  ideal  gate  of  uniform  sheet  resistance,  the  variation  of 
gate  potential  distribution  is  shown  in  Fig.  4.16  (The  simulation  is  not  self  consistent, 
merely  to  indicate  the  effect  of  oscillation  on  the  potential  along  the  gate).  The  larger  the 
sheet  resistance,  the  larger  the  variation  of  voltage  along  on  the  gate. 

The  upper  limit  of  sheet  resistance  Rs  can  be  estimated  by  setting 


(4.7) 

where  Iq  is  the  average  gate  current,  AIq  is  the  peak-to-pcalc  variation  of  Ig.  Wq  is  the  gate 
width,  and  is  applied  electric  field.  That  way,  the  channel  potential  will  retain  its 
linearity  and  the  electric  field  can  be  kept  constant  to  ensure  that  multiple  domains  occur.  In 
the  absence  of  diode  current  (no  spillover),  the  variation  of  gate  current  is  given  by  the 
variation  in  capacitor  current  ic: 


Therefore, 


A,  •  UUA 

AIq  -  1q  -  qWAn,Vj  «  — ^ 

*'s 


Here,  vj  *Ay  /At  is  the  domain  velocity,  W  is  the  channel  width,  and  dn,  is  the  peak-to- 
peak  sheet  electron  density  in  the  channel.  The  final  result  for  the  upper  lound  of  the  gate 
sheet  resistance  is; 

— <-3r> 

(4.10) 

Taking  data  from  the  computer  simulation  result  for  a  device  having  channel  doping 
ND»I.2xlO*^cm'^,  Anpp«7xl0ll  cm'2,  vd«10^cm/s,  the  upper  limit  for  Rs  is  equal  to  7.1 
kfl  per  square  for  an  electric  field  of  80(X)V/cm. 

Too  low  a  resistivity  results  in  an  unacceptable  amount  of  heating  in  the  gate  and  higher 
resistance  allows  the  gate  voltage  to  be  moduled  locally  by  electron  domains.  To  be  on  the 
safe  side,  Rs»3kn  per  square  will  be  the  target  for  choosing  proper  gate  material. 
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Therefore,  the  present  existing  RG-MESFET  oscillators  need  to  be  modified  to  have  a 
lower  gate  resistance  to  be  able  to  screen  the  space  charge  in  the  channel. 


Figure  4. 16  Gate  voltage  variation  due  to  the  existence  of  charge  packets  in  the  channel 
of  an  RG-MESFET  with  gate  resistance  of  (a)  30kfl/square,  and  (b)  3 
kfl/square  [24]. 
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4.5  CrSiO  Film  as  Resistive  Gate  Material 

As  discussed  in  the  previous  section,  the  single  domain  generation  and  propagation  in 
the  RG-MESFET  oscillator  is  caused  by  the  high  resistivity  of  the  gate.  In  order  for  the 
contiguous  domain  oscillation  mode  to  occur,  the  sheet  resistance  of  the  gate  has  to  be 
around  3kf2/square.  The  present  resistive  gate  is  made  of  WSiN  and  has  a  thickness  of 
1  lOOA  and  a  sheet  resistance  of  30kG/square. 

To  modify  the  existing  RG-MESFET  oscillator,  it  is  necessary  to  etch  away  the  original 
resistive  gate  and  deposit  a  new  resistive  gate  with  the  right  sheet  resistance.  Keep  in  mind 
that  since  we  also  want  to  eliminate  the  gate  overlapping  on  the  source  and  drain,  new  gate 
contacts  need  to  be  formed  along  with  a  new  gate  film.  There  are  two  ways  to  accomplish 
this.  Fig.  4.17  (a)  shows  the  sequence  of  depositing  the  new  gate  first  and  then  depositing 
the  contact  metal.  This  process  requires  two  lithography  steps.  Fig.  4.17(b)  is  just  the 
opposite,  but  requires  three  lithography  steps. 

The  new  gate  film  has  to  make  a  low-leakage  Schottky  contact  to  the  semiconductor.  It 
has  to  be  very  uniform  to  provide  a  uniform  electric  field  in  the  RG-MESFET  channel.  The 
thickness  should  be  from  1000  to  2000A.  With  a  resistivity  one  order  of  magnitude  smaller 
than  the  original  WSiN  gate,  the  sheet  resistance  of  the  material  should  be  from  3kQ/square 
to  1.5kG[/square.  In  addidon,  an  ohmic  contact  is  required  for  this  film  with  a  very  small 
contact  resistance. 

Song  et  al.[26]  have  fabricated  an  electron-beam  evaporated  Cr-SiO  /GaAs  diode  and 
measured  the  current-voltage  characteristics.  According  to  their  result,  a  very  good 
Schotticy-diode  characteristics  for  the  Cr-SiO  cermet/GaAs  diode  can  be  obtained  by 
annealing  at  425®C.  The  leakage  current  is  nearly  two  orden  of  magnitude  lower  than  those 
of  annealed  A1  and  Cr  Schottky  diodes.  Moreover,  the  annealing  temperature  and  the 
composition  of  the  film  were  found  to  be  useful  parameters  for  adjusting  the  values  of  the 
resistivity  [27-30].  Therefore,  Cr-SiO  cermet  film  is  the  choice  for  noodifying  the  RG- 
MOSFET. 

As  indicated  by  Song  et  al  [26],  the  evaporation  temperature  of  Cr  and  SiO  are 
relatively  low  and  close  to  each  other,  which  makes  the  Cr-SiO  cermet  a  good  choice  for 
electron-beam  evaporation.  The  use  of  a  mixture  of  powders  as  an  evaporation  source 
makes  it  convenient  to  change  the  film  composition. 

However,  due  to  the  different  fabrication  conditions  used  and  the  varying  nature  of  the 
cermet  film,  the  resistivity  of  the  resulting  films  needs  to  be  carefully  controlled  for  this 
panicular  application.  There  has  been  no  report  on  the  fabrication  of  the  ohmic  contact  on 


New  gate  contacts 


Figure  4.17  Modification  of  the  RG-MESPTET  oscillator. 

The  gate  of  Fig.  4.1  has  to  be  etched  away  first.  Two  process  sequences 
can  ^  used,  (a)  deposit  the  new  gate  first  and  then  the  gate  contacts,  (b) 
deposit  the  contact  metal  first  and  then  the  new  gate  film. 


the  cermet  films,  which  is  very  critical  to  us  and  many  other  applications,  especially  when 
the  sheet  resistance  is  low. 

To  investigate  the  resistivity,  60  wt.%Cr-40  wt.%SiO  and  55  wt.%Cr-45  wt.%SiO 
resistors  were  fabricated.  Heat  treatment  was  performed  the  same  way  as  it  would  be  in  the 
modification  of  the  RG-MESFET  oscillators  for  the  60  wt.%Cr-40  wt.%SiO  cermet  film. 
Different  metalization  schemes  were  investigated  in  order  to  reduce  the  contact  resistance. 
To  better  control  the  sheet  resistance  of  the  film,  plasma  etching  using  CF4  and  O2  was 
used.  The  RG-MESFET  oscillators  were  modified  by  etching  the  original  gate  away. 
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depositing  Cr-SiO  cennct  film,  and  contacting  it  with  Ti/Au  metal  layers.  IV  measurements 
indicate  that  the  cermet  films  basically  meet  the  requirements  stated  above. 

4.5.1  Fabrication  of  Cermet  Resistors 

The  source  of  the  evaporation  was  prepared  by  mixing  powders  of  Cr  (50  mesh)  and 
SiO  (120  mesh)  (both  were  from  Johnson  Matthey  Inc,  Alpha  Products,  152  Andover 
Street,  P.O.  Box  299,  Danvers,  MA  01923)  into  60  wt.%Cr-40  wt.%SiO  and  55  wt.%Cr- 
45  wt.%SiO  compositions. 

Si  substrates  with  thermally  grown  oxide  films  of  about  1500A  thickness  (not  critical, 
mainly  for  insulating  purposes)  were  prepared  by  degreasing  in  trichloroethylene,  acetone 
and  methanal  with  ultrasonic  agitation  for  5  mins  each  and  rinsing  with  DI  water.  Positive 
photoresist  (AZ1350J-SF)  was  spun  at  40(X)  r.p.m.  for  30  seconds  and  optically  patterned 
in  preparation  for  lift-off.  The  pattern  includes  resistors  with  different  lengths  and  widths 
as  well  as  1cm  x  1cm  squares.  The  former  are  for  measuring  the  film  sheet  resistance  and 
contact  resistance  by  probing  the  contact  metal.  The  latter  are  for  measuring  the  sheet 
resistance  measurement  by  4-point  probe.  The  resistors  have  two  metal  pads  (200)im  by 
200|im)  on  the  ends,  also  formed  by  lift-off  process. 

The  evaporation  of  the  cermet  film  is  carried  out  in  a  Varian  E-beam  evaporation  system 
with  a  deposition  rate  and  thickness  monitor.  The  beam  voltage  of  the  electron  gun  is  about 
6  kV,  which  is  fixed  for  our  system.  Because  of  the  high  electron  gun  voltage,  the  beam 
current  is  very  small,  around  0.01  A.  The  resulting  deposition  rate  is  around  20-40  A/s.  The 
system  pressure  before  evaporation  is  typically  2x10*'^  torr.  The  evaporation  is  monitored 
to  stop  when  the  sheet  resistance  reaches  the  desired  value. 


4.5.2  Film  Characterization 

Film  thickness  is  measured  by  an  Alpha-Step  200  manufactured  by  Tcncor 
Instruments,  which  is  a  computer  controlled  stylus  instrument  Measurement  is  carried  out 
by  physical  contact  of  stylus  and  films. 

The  sheet  resistance  is  measured  in  two  ways.  The  easy  and  fast  one  is  by  a  Digital 
Resistivity  Test  System  (4-poir.t  probe).  The  other  way  is  to  probe  the  resistors  with 
different  length  and  construct  an  RW  vs.  L  plot  Here  R  is  the  measured  resistance  of  a 
resistor  with  length  L  and  width  W.  The  slope  of  the  straight  line  is  the  sheet  resistance  of 
the  film  and  the  intercept  derived  by  W  is  the  total  contact  resistance  2Rc,  since 


R-R,^  +  2R, 
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(4.11) 


4.5.3  Result  and  Discussion 


4.5.3. 1  Sheet  Resistance  of  Cermet 

The  composition,  thickness,  sheet  resistance,  and  bulk  resistivity  of  the  films  at  room 
temper  before  the  thermal  annealing  are  shown  in  Table  4.3. 


Table  4.3  Properties  of  the  ev  api^ratcd  cermet  films 


Cr  wt.  % 

Sample 

Thickness 

(X) 

R,  at  30O^K 
(kfi'squarc) 

P  (xlO-3) 
(Q-cm) 

55  % 

CR3 

3983 

38.88 

1.55 

CR9 

3433 

54.61 

1.87 

60% 

CR13 

2523 

9.0-4 

0.23 

CR14 

2230 

5.57 

0.12 

CR15 

3440 

6.42 

0.22 

CR16 

1513 

5.93 

0.09 

BlOF 

2313 

5.22 

0.12 

Note  that  the  data  arc  scattered.  This  is  because  the  source  is  made  of  powder  material 
with  high  melting  points.  The  source  does  not  melt,  only  sublimes  during  the  evaporation 
[32].  One  component  might  evaporate  a  little  faster  than  the  other.  Therefore,  during  the 
evaporation,  it  is  important  that  the  timing  of  each  action  (mm  on  high  voltage,  turn  on 
emission  current,  open  the  shatter,  etc.)  is  kept  the  same  for  every  run  to  ensure 
reproducibility.  Also,  fresh  material  is  desired  for  each  evaporation.  One  better  way  is  to 
compress  the  powder  material  so  that  the  particles  arc  in  contact  to  each  other  and  the  heat 
transfer  is  more  uniform.  As  seen  in  the  table,  the  resistivity  of  60  wt.%Cr-40  wt.%SiO 
cermet  is  a  better  choice  for  the  RG-MESFET  device  modification.  Although  the  resistivity 
is  bigger  than  the  desired  value  of  0.03  O-cm,  it  is  shown  later  that  thermal  a.’'nealing  will 
reduce  the  resistivity  considerably.  A  thickness  of  around  2000A  is  desirable  since  the 
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relative  thickness  variation  is  smaller,  which  is  important  for  the  CDO  deWcc,  The  surface 
profile  of  the  resistor  from  the  Tencor  Alpha-step  measurement  indicates  that  the  film  is 
fairly  uniform. 

Annealing  of  the  cermet  films  is  performed  in  a  Marshal  furnace.  The  experimental 
results  are  summarized  in  Fig.  4.18.  When  exposed  to  the  air,  the  sheet  resistance  is  higher 
compared  with  that  under  otherwise  the  same  annealing  conditions  due  to  the  existence  of 
oxygen.  In  all  cases,  reduced  sheet  resistances  arc  measured  after  annealing.  Continued 
annealing  will  further  reduce  the  resistivity  of  the  film. 


Annealing  Time  (mins) 


Figure  4.18  Relative  resistance  of  Cr-SiO  films  (60  wt.%Cr-40  wt.%SiO)  as  a  function 
of  annealing  time.  The  sheet  resistance  is  measured  by  4-point  probe. 


4.5.3.2  The  Contact  Resistance 

The  contact  resistance  between  the  Cr-SiO  film  and  the  metal  pad  is  an  important 
parameter  that  has  to  be  kept  small. 

Experimental  results  show  that  a  plasma  clean  is  necessary  right  before  the  cermet 
evaporation.  This  is  seen  in  Fig  4.19  for  resistors  with  Ti/Au  metal  pads  (500  A/IOOOA, 
evaporated  in  Varian  E-beam  evaporation  system).  The  Ti/Au  pads  were  evaporated  first.  If 
the  preclean  is  omitted,  the  resistance  of  die  resistor  varies  considerably  ( Fig.  4.19  (a)). 
Although  the  average  slope  is  close  to  the  sheet  resistance  from  the  4-point  probe 
measurement,  the  contact  resistance  is  fairly  large  and  nonuniform  from  resistor  to  resistor. 
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hnlran  makri  the  contact  trtiiitanre  amailer  and  more  untfom  (f  ig.  4,19  (b)).  It  n 
iKbrvfd  to  irnKHT  the  rcstdiial  resiat  at  (he  surface  of  the  tnctal  pad.t.  Ilte  ptrclcan  is  done 
u^iiig  O  4  in  a  Tci  hnics  I’l!  IIA  plastna  system  for  10  sec.  with  RF  power  of  40  W  and 
chamber  pressure  of  4(X)  inTorr. 


0  200  400  floo  Sty)  1000  i:oo 

L(kim) 


I  igtire  4  19 


400  600  800  KXy)  1200 

L(kun) 


I^ffei  t  of  preclean  on  the  contact  resistance.  Metal  is  cvafiorated  beF  re  the 
cermet  film  Fach  sytrihol  represents  a  set  of  rcsiston  of  different  length  in  a 
die  The  solid  straighi  line  is  ihe  rcsuli  measured  fmm  the  4  point  prof>e 
(a)  No  prrcleanmg 

(h)  IhTclean  using  (T4  in  »  Technics  plasma  system  with  RF  power  of  40 
W  amj  chamber  pressure  of  4(X)  mTorr. 
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As  seen,  prrclcan  greatly  impmves  the  quality  of  the  contacts. 

If  the  order  of  metal  pad  and  certnet  film  is  interchanged,  i.c.,  Ti/Au  is  placed  on 
cennet,  k)\ver  contact  resistance  is  observed  (Mg.  4  20).  litis  is  because  when  the  cermet 
is  place  on  Ti/Au  pul,  the  inert  nature  of  the  gold  prevents  adhesion  to  the  cemtet  layer  by 
chemical  Ixitiding.  The  transition  metal  titanium  promotes  the  adhesion  to  both  ccmx’t  and 
gold. 


L  (Jim) 

Figure  4.20  Rffcct  of  the  fabrication  sequence  of  making  resistors.  Each  symbol 
represents  a  set  of  resistors  of  different  length  in  a  die.  The  solid  straight 
line  is  the  result  measured  from  tltc  4*point  probe. 

(a)  Cermet  on  Ti/Au. 

(b)  Ti/Au  on  cermet. 
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Case  (b)  has  a  smaller  contact  resistance  and  less  scatteied  than  that  of  case 
(a). 

In  order  to  minimize  the  contact  resistance,  several  metal  interconnections  were 
investigated.  All  the  resistors  are  fabricated  by  evaporating  metal  on  top  of  the  cermet. 
Prcclean  is  performed  before  the  metalization.  Ti,'Au  and  Ni/Ti/Au  arc  evaporated  in  a 
Varian  E-bcam  system  'vith  thickness  monitor.  A1  and  Cr/Au  arc  evaporated  in  an  NRC 
system.  Thickness  is  not  well  controlled.  Fig.  4.21  shows  the  effect  of  four  types  of 
metalization  on  the  contact  resistance.  Trying  to  obtain  a  better  contact,  a  4(X)®C  annealing 
in  N2  for  10  minutes  is  performed  on  all  the  resistors.  As  shown,  better  ohmic  contacts  are 
obtained  after  annealing  except  for  Cr/Au  metalization. 


LOim) 


Figure  4.21  Effect  of  Al,  Cr/Al,  Ni/Ti/Au,  and  Ti/Au  metalization  on  the  contact 
resistance.  The  resistance  of  all  the  samples  is  measured  before  (open 
symbols  and  dased  lines)  and  after  (solid  symbols  and  solic  lines)  thermal 
annealing  at  400^0  in  N2  for  10  min. 
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After  annealing,  the  contact  resistances  are  all  reduced.  The  data  arc  less 
scattered  after  annealing  except  for  the  Cr/Au  meulization.  For  Ti/Au 
metalization,  the  contact  resistance  is  eliminated  entirely. 
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Figure  4.21,  continued. 


4.5.3.3  Plasma  Etch  of  the  Cermet  Film 

It  is  desirable  to  find  an  etchant  for  the  cermet  film  if  necessary.  The  etchant  gas  of 
CF4, 02,  N2  and  the  mixture  of  CF4  and  O2  were  tested.  Only  the  mixture  of  CF4  and  O2 
is  able  to  etch  the  cermet  film.  The  etch  rate  is  453A/min  using  IOOCF4  :  5OO2  at  a  chamber 
pressure  of  400  mTorr  and  power  of  150W  in  a  Branson/IPC  barrel  etcher.  In  contrast,  the 
etch  rate  of  Au  under  the  same  conditions  is  80A/min. 


4.5.4  RG-MESFET  Device  with  60%Cr-40%SiO  Cermet  Gate 

The  result  of  CrSiO  film  fabrication  indicate  that  the  modification  procedure  depicted  in 
Fig.  4.17  (a)  is  a  good  choice.  It  requires  only  two  lithogartphy  steps.  With  the  Ti/Au  metal 
on  top  01  the  CrSiO  gate,  the  contact  resistance  is  minimized. 

60  wt.%Cr-40  wt.%SiO  cermet  film  was  used  to  form  a  resistive  gate  on  an  RG- 
MESFET  oscillator.  To  do  that,  PMMA  positive  photoresist  was  spun  on  the  original  RG- 
MESFET  device  at  400  r.p.m.  for  30  sec.  and  baked  at  160®C  for  4  hours.  A  window 
exposing  old  WSiN  was  defined  by  electron-beam  direct  writing.  WSiN  was  etched  away 
using  CF4  in  a  Branson/IPC  Barrel  etcher  with  power  of  150W  and  chamber  pressure  of 
0.4  torr.  A  60  wL%Cr-40  wt.%SiO  cermet  film  was  then  evaporated  and  lifted  off.  The 
second  lithography  defines  the  metal  contacts  and  Ti/Au  metalization  to  the  cermet  was  e- 
beam  evaporated.  As  a  result  of  the  4-hour  bake  at  160®C  and  the  heat  exposure  during  the 
metalization  in  this  sequence,  the  resistivity  of  the  resulting  cem^t  was  stabilized.  After  the 
metalization,  a  10  min  anneal  at  400^C  was  done  to  reduce  the  contact  resistance.  The 
resulting  cermet  gate  had  a  sheet  resistance  of  0.5  kG/square.  The  TV  curve  is  shown  in 
Fig.  4.22I  The  fact  that  the  gate  current  Ig2  varies  linearly  with  respect  to  the  applied  gate 
voltage  Vgg  is  a  strong  indication  that  the  leakage  current  of  the  ceimct/GaAs  Schottky 
junction  ii  small. 

Unforttinately,  because  the  gate  resistance  was  so  small,  the  gate  was  destroyed  before 
any  oscillation  was  observed.  Since  this  modification  process  has  shown  that  CrSiO  can 
reduce  the  gate  resistance  and  form  low  leakage  Schottky  contact  with  GaAs,  it  is  ideal  to 
design  new^G-MESFET  device  and  incorporate  the  CrSiO  resistive  film  process  in  the 
fabrication  skjuence. 


Gate  Voltage  Vqq  (V) 


Figure  4.22  IV  characteristics  of  a  RG-MESFET  device  with  a  60  wt.%Cr-40  wt.%SiO 
cennet  gate.  LxW=l(hc20. 


4.6  Conclusion 

The  resistive  gate  MESFET  devices  have  been  fabricated  to  test  the  contiguous  domain 
concept  The  device  with  channel  doping  density  of  1.2x10^’^  cm'3  and  channel  depth  of 
ISOOA  exibits  microwave  oscillation.  The  oscillation  is  identified  as  single  domain 
oscillation  due  to  the  high  gate  sheet  resistance  of  SOkH/square.  Preliminary  analysis 
indicates  that  the  gate  sheet  resistance  has  to  be  lowed  to  about  3kG/square  to  avoid  the  gate 
potential  variation  caused  by  the  channel  oscillation.  Cermet  fabrication  process  was 
developed  to  obtain  a  desirable  resistive  film  for  future  fabrication  of  CDO  devices. 


CHAPTER  5 

EQUIVALENT  CIRCUIT  SIMULATION  OF  THE  RG-FET 


5.1  Introduction 

As  shown  in  the  last  chapter,  single  domain  oscillation  is  obtained  in  the  RG-MESFET 
device.  However,  the  TOF  program  is  not  able  to  simulate  the  single  domain  since  in  the 
TOF  program  the  channel  potential  is  only  determined  by  the  gate  voltage.  The  effect  of 
source  and  drain  is  neglected.  In  order  to  simulate  contiguous  domain  and  single  domain 
oscillation,  a  very  simple  lumped  circuit  model  is  developed.  This  model  demonstrates,  for 
the  first  time,  the  oscillation  mode  transition  in  a  doped  channel  RG-FET  structure. 

The  equivalent  circuit  model  is  established  through  the  equations  of  state  approach  [33]. 
In  this  approach,  one  first  derives  the  current-voltage  relationships  (or  in  some  instance, 
charge-voltage  relations)  which  describe  the  internal  operation  of  the  device.  Next  an 
equivalent  circuit  is  deduced  by  combining  idealized  circuit  elements  in  such  a  manner  as  to 
yield  the  same  current- voltage  (or  charge-voltage)  equations. 

The  equations  of  state  include  throe  basic  relationships  expressing  continuity  in 
semiconductor  devices.  The  change  ir  hole  density  in  a  region  over  a  interval  of  time  must 
equal  the  excess  of  flow  into  the  region  for  that  interval  plus  the  excess  generation  of  hole- 
electron  pairs  during  the  interval.  A  similar  relationship  applies  to  electron  density.  In 
addition,  the  changes  in  potential  of  a  region  are  related  to  the  charge  through  Poisson's 
equation. 

Assuming  no  generation-recombination  sources  inside  the  semiconductor,  the  three 
continuity  equatii  sare: 


9Py 

dt 

9nv 

"dr 


-  e,V^(l>  =  q  ( Pv  -  n,  +  Njj  -  N^  ) 


(5.1a) 

(5.1b) 

(5.1c) 
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where  pv  is  the  volume  density  of  holes,  ny  is  the  volume  density  of  electrons,  and  ^  is  the 
intrinsic  Fermi  potential  in  the  channel.  From  Boltzmann  transport  theory,  the  total  electron 
and  hole  current  density  (AC  +  DC),  Jn  and  jp  in  equation  (5.1)  can  be  wiitt<*n  as  functions 
of  quasi-Fermi  levels  for  electrons  and  holes,  vn  and  vp: 


jp  =  -  qiipPvVvp 

jj^  =  -q|i„nvVvj^ 


(5.2a) 

(5.2b) 


or  alternatively,  Jn  and  jp  can  be  written  as  functions  of  ()>,  ny  and  py,  consisting  of  drift  and 
diffusion  component: 

jp  =  Pv  -  qDp  ^Pv  (5.3a) 

jN  =  ‘l^^l^"v  +  qD„VI^ 


where  )Xn  and  4p  are  the  electron  and  hole  mobilities  and  Dn  and  Dp  are  the  electron  and 
hole  diffusivities.  4pt  and  Dp  are  functions  of  electric  field 

From  Eqn  (5.1),  we  can  derive  an  equivalent  circuit  model  for  specific  device 
structures.  We  will  work  our  way  through  from  a  very  simple  one-dimensional 
semiconductor  stub  to  a  doped  channel  RG-FET  struemre  to  show  the  simplicity  and 
validity  of  the  model 


5.2  Equivalent  Circuit  Model  for  a  Semiconducter  Stub 

A  semiconductor  stub  is  chosen  to  illustrate  the  ideas  regarding  lumped  circuit  model. 
The  structure  to  be  simulated  is  shown  in  Fig  5.1(a).  The  problem  is  one-dimensional, 
assuming  all  the  quantities  are  uniform  along  the  x  and  z  directions.  Equations  (5.1)  and 
(5.2)  then  become: 


where 


,  1  ^jp 

dt  q  ay 


=  0 


at  'q  dy 
d^A 

-  ^  Pv  -  "v  +  Np  -  ) 


jfr  *  qPn(  *  ^  ^ 


(5.4a) 

(5.4b) 

(5.4c) 


(5.5a) 


(5.5b) 
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Choosing  the  potential  reference  to  be  the  Fermi  level  at  the  ca±ode  (Fig.  5.1b),  we  define: 


FJcathode)  -  FpfjO 

Vp(y)=-Jl - - - 1— 

~  q 

Fp,(cathode)  -  E.(y) 
4,(y)  =  -S - - - ! - 


(5.6a) 

(5.6b) 

(5.6c) 


Hole  and  electron  volume  density  can  be  expressed  in  terms  of  vp,  vn  and  <j>  using 
Boltzmann  statistics, 

(E-FJ/kgT  q(Vp-«)/kBT 
Pv*nje  =nje 

The  time  variation  of  pv  and  ny  can  be  found  by  differentiating  (5.7): 

dn 


_  Q  ^  d  ( jf.  „  \ 

Dt  "k^T  "dt^^  N> 


Substituting  (5.8a  &  b)  in  (5.4a  &  b),  we  have 

i^P.-|(Vp-«+^  =  0 

5^n.|'(v„-«  +  ^=0 

The  time  derivative  of  (5.4c)  yields: 

.2 


a  d  <|»  „  .  3py  9nv 
■  .  •»  “  4  t.  "3:  57" 


dy 


dt  dt 


After  defining  a  displacement  current 


_  d<}> 

■  Ay  dt 


(5.10)  becomes: 


(5.7a) 

(5.7b) 


(5.8a) 

(5.8b) 


(5.9a) 

(5.9b) 


(5.10) 


(5.11) 


(5.9c) 
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The  govcrnir  ferential  equations  (5.9a,  b  &c)  result  in  a  disaibuted  model  in 
nature.  The  din  s  of  the  element  are  infinitesimal.  The  process  of  replacing 
differential  equations  by  difference  equations  turns  the  infinitesimal  elements  into  finite 
lumps.  In  this  case,  the  semiconductor  stub  is  divided  into  m  small  sections.  For  each 
section  of  lengtli  Ay,  a  single  value  of  the  pertinent  variable  -  hole  density  pv,  electron 
density  ny,  and  static  potential  »i»  -  is  selected  to  represent  the  whole  element  Equations 
(5.9)  become: 


Define 


2 

■j^  Pv  Ay  ^  (Vp-  9 )  =  jp(y  -  Ay)  -  jp(y) 

^  n.  Ay  ^  (v^-  9 )  =  jN(y  -  Ay )  -  j^Cy) 
B 


(5.12a) 

(5.12b) 


_ql 


2 

Pv  Ay  ^(  9  -  Vp)  +-j^ny  Ay  ^(  9  -  Vj^)  =  ji(y  -  Ay)  -  jj^(y) 

^  (5.12c) 


T^Pv  Ay 

(farads/cm^) 

(5.13a) 

q^  A 

e. 

(farads/cm^) 

(5.13b) 

> 

Ay 

(farads/cm2) 

(5.13c) 

The  equations  of  the  state  for  the  semiconductor  stub  become 


Cp  ^  ( Vp-  9 )  =  jp(y  -  Ay )  -  jp(y) 

c^|-(Vf^-9)=jN(y-Ay)-jN(y) 

CV^^(9-  Vf,)  +  Cp^(9-  Vp)  =jj(y  -  Ay  )  -jj(y) 


(5.14a) 

(5.14b) 

(5.14c) 


The  form  of  Equation  (5.14)  clearly  suggests  node  equations.  It  can  be  easily 
interpreted  in  terms  of  the  three-line  equivalent  circuit  model  (Fig.  5.1c).  The  top  and 
bottom  lines  can  be  thought  of  as  representing  the  conduction  and  valence  band, 
respectively.  The  central  line  represents  the  displacement  current  path  [33]. 

Although  Fig.  5.  Ic  looks  like  a  small  signal  equivalent  circuit,  the  definition  of  Cp  and 
Cn  are  in  terms  of  the  total  hole  and  electron  densities.  Using  this  circuit  correctly,  we  can 
simulate  the  large  signal  response. 


tJiJJI  .^WUIfipwipiBpWIIIWyfWW^^  HIJti||!M|;nii||ii  II  « 

lamai  ^  ^'>ia,.'sai.»,.i., 
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5.3  Equivalent  Circuit  Mtxlc  1  for  a  Gunn  DukIc 

5.3.1  Simulation  I*r(X'C(1ure 

In  this  section,  a  procedure  for  simukting  a  Gunn  Dicxle  (an  n  typ«  semiconductor 
stub)  using  the  established  equivalent  circuit  is  described,  which  is  applicable  at  a  later 
point  in  the  simulation  of  a  doped  channel  RG-MiT.  Since  we  arc  only  interested  in  the  n 
type  semiconductor,  the  p-line  is  dntpped  out  of  the  pic  ture.  The  corresponding  equivalt  nt 
circuit  is  shown  in  Fig.  5.2. 


Given  tSie  instantaneous  electron  distribution  and  the  terminal  l>oundary  conditions,  the 
channel  potentitl  $  can  be  found  by  solving  a  one -dimensional  Poisson's  equation.  In 
order  to  get  Poisson's  equation  from  the  equivalent  circuit,  we  need  to  work  with  the 
charge  balance  at  each  node  in  the  Mine.  At  node  j,  the  uwal  charge  on  all  capacitor  plates  is 
zero. 


0 


(5.15) 


Notice  that  all  the  charges  Q  in  the  above  equation  are  large  signal  quantities.  Since  Gg  is  a 
constant  capacitor,  the  total  charge  it  simply  the  prcxluci  of  the  capacitance,  the  voltage 
across  it,  and  the  related  area; 


IM 


1  .i,iHi[ii>^m'iwyyii^sFy.'  >in.iiii,.ii.i,.i.i.i.iijijijiii 
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uhcrr  W  and  M  are  the  width  and  height  of  the  stub,  rejpcctively.  For  the  voltage 
dcf'endcnt  Cn.  the  total  charge  is  obuined  by  integrating  tltc  differential  charge  Cn  dfO-vs,) 
over  the  voltage  ddTcrrnce; 


(♦■»»,) 


Q^„wn  j  ^-Ayn,d(0  Vj^) 
’n'o  ® 

-Wn3:jrAyn.  J  e 


»  WII  q  Ay  n,  (  e 


d((>-v^) 

qt^  v^iAjT  q<* 


) 


(5.17) 


where  (0  vn)o  >*  determined  by  the  DC  elect  jii  charge  on  the  capacitor  through  the 
rclatuanship 


Njj-n,e 


<»(♦•  Vfl'SiT 


Therefore, 


ONj-WlIqAy  (n^^-  Nj>) 

Eeptation  (5.15)  now  becomes: 

C^(0..,-20,>0j^,)-q(n,j-  Np)  Ay 
Recalling  that  Ck  "  fi  /  Ay.  the  final  charge  balance  at  ncxic  j  is: 


(5.18) 

(5.19) 

(5.20) 


Ay' 


'  q  (n,j-  Np) 


(5.21) 


This  is  exactly  the  same  as  the  discretized  Poisson’s  equation 

d^6  /  vr  \ 

e,-^-q(n,-Np) 

‘ly  (5.22) 

To  simulate  the  Gunn  oscillation,  the  device  is  divided  into  m  sections.  (5.21)  has  to  be 
solved  numerically  for  the  channel  potential  0  with  the  appropriate  boundary  conditions, 
given  the  electron  density  ny.  If  ohmic  contacts  arc  assumed  for  the  sourcc/drain,  the 
excess  charges  on  the  two  contacts  are  zero.  Tlicrcforc.  the  volume  densities  at  the  source 
and  drain  are  fixed  at  the  channel  doping  density,  and  so  arc  the  potentials  at  the  contacts. 
The  boundary  conditions  for  Poisson's  equation  arc: 
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^  In  (-J1) 

knT  Nn 


(5.23a) 


(5.23b) 

where  nj  is  the  intrinsic  density  of  the  semiconductor.  Once  the  channel  potential  is 
obtained,  the  electric  field  in  the  chanrxl  is  determined  by 


^  (5.24) 

The  nonlinear  velocity- field  curve  of  (3.10c)  is  used  to  determine  the  velocity  of  electrons 
at  each  grid  point: 


'‘■I-- 


^1‘c'V’' 


-.1/2 


Mo'V 


(5.25a) 


■  Ay  e 


1  + 


(5.25b) 


and  the  electron  current  flow  j'n  can  be  calculated  from  Eqn  (5.5b).  From  the  current 
distribution,  the  time  derivative  of  the  electron  concentration  is  obtained  from  the  continuity 
equation  (5.4b).  From  this  time  derivative,  the  electron  density  an  instant  later  is  computed: 


with  the  boundary  conditions: 


(5.26) 

(5.27a) 

(5.27b) 


and  the  cycle  repeats  [34).  The  same  algorithm  shown  in  Fig.  3.3  is  used  here  as  in  the 
Time-Of-Flight  program.  The  only  difference  is  in  the  calculation  of  the  potential 
distribution.  Here  a  set  of  linear  equations  has  to  be  solved  simultaneously. 
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5.3.2  Simulation  Results 

Fig.  5.3  shows  the  simulation  results  for  a  Gunn  diode  biased  at  2V,  using  the  above 
formuation.  The  device  has  a  length  of  2^m,  width  of  l|im  and  height  of  lp.m^  The  doping 
density  is  1x10*^  cm'3.  The  number  of  grid  points  is  256.  The  typical  time  step  At  is  in  the 
femto  second  range  and  is  self-adjusting.  The  nonlinear  velocity  field  curve  is  used  to 
obtain  the  oscillation.  The  parameters  in  (5.25b)  arc:  )ao*52(X)  cm^/V-s,  Av=9. 55x10'^, 
Bv»  1 .3 1 6x  107,  Cv«5 .8x  lO’,  Dv«6.788x  105,  and  Fv=0.3 1 6. 

At  time  t=0,  the  surface  electron  density  (ns=nvH  at  the  value  NqH)  is  constant  except 
in  the  middle  of  the  device  where  an  arbitrary  dipole  domain  is  intJ-oduced  (Fig.  5.3(a), 
black  dots).  As  time  evolves,  this  dipole  grows,  expands  and  quickly  stabilizes.  This 
stabilized  dipole  travels  with  a  velocity  of  9.15x10*^  cm/s  toward  the  drain,  as  determined 
from  Fig.  5.3(a). 

The  applied  electric  field  is  V(yL=10kV/cm  without  the  dipole.  With  the  existance  of  the 
introduced  small  dipole,  the  electric  field  inside  the  dipole  increases  and  the  field  outside 
decreases  a  little.  As  the  dipole  grows,  the  field  inside  enhances  and  the  field  outside 
decreases  further,  as  shown  in  Fig.  5.3(b).  The  same  effect  can  also  be  seen  in  Fig.  5.3(c), 
where  the  slope  of  the  curves  represents  the  electric  field. 


I 

I 


♦  (V)  ^(kV/cm) 
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5.4  Equivalent  Circuit  Model  for  a  Doped  Channel  RG-FET 

The  cross  section  of  a  doped  channel  RG-FET  in  Fig,  5.4  shows  the  device  to  consist 
of  a  p-type  substrate  with  an  n-diffusion  overlaid  by  an  insulator  and  a  resistive  gate.  The 
n-diffusion,  referred  to  as  the  channel  region,  looks  like  a  "semiconductor  stub"  with  the 
assumption  that  it  is  uniformly  doped. 


Figure  5.4  Cross  section  of  a  doped  channel  RG-FET. 


5.4.1  Effect  of  the  Gate  on  the  Equivalent  Circuit  Model 

Normally,  the  substrate  is  very  lighdy  doped,  and  the  depletion  region  in  the  channel 
due  to  the  channel-substrate  junction  can  be  treated  as  zero.  In  the  case  of  the  gradual 
channel  approximation  [13],  the  cross  section  in  the  x  direction  should  look,  in  general, 
like  a  regular  doped  channel  MOSFET  (Fig5.5a).  Current  flow  in  the  x  direction  is  zero. 


The  only  governing  equation  is  Poisson's  equation  in  x  direction: 


q  ( "v  -  NjP 


(5.28) 


Fig.  5.5(c)  is  the  band  diagram  obtained  by  solving  (5.28)  for  the  block  charge  diagram 


shown  in  Fig.  5.5(b).  The  channel  is  connected  to  the  gate  through  two  capacitors,  Cox  and 


Cd-  Cd  is  the  depletion  capacitance  associated  with  the  depletion  region  of  width  x^j.  Cox  is 


the  insulator  capacitance.  The  voltage  across  Cd  can  be  determined  by  the  following 
relationship: 


<^+Ayi 


Cox  Qj 


Figure  5.5 


(a)  Cross  section,  (b)  block  charge  diagram,  (c)  band  diagram,  and  (d) 
equivalent  circuit  model  of  a  dopcd-channel  MOSFET  structure  along 
die  X  direction,  under  the  gradual  channel  approximation. 


96 


f^-AV  --I-  -AV„+^ 

>q  d  q  ox  q 


-  -  =  0 
1  o 


2q  q  q  o  (5  29) 

where  <l>'ni  and  x'  are  the  effective  metal  workfunction  and  electron  affinity  of  the 
semiconductv  ,  both  relative  to  the  conduction  band  of  the  oxide.  Assuming  that  the  gate 
charge  is  equal  to  the  semiconductor  bulk  charge,  the  voltage  drop  on  the  oxide  is: 


where 


Define: 


V 

'-ox 


V’  -Y*-— 

Substituting  (5.30)  and  (5.32)  into  (5.29),  we  have: 

•‘Vj  +  yiVjAVj  -('1>  +  v;-Vq)  =  0 
Solving  (5.33),  the  potential  associated  with  the  band  bending  AVj  is: 


(5.30) 


(5.31) 


(5.32) 


(5.33) 


AV,.4,  +  v;.V<j,-V„.7vJ-f2V,(4.*v;.VjP 


(5.34) 


From  Fig.  5.5(d),  it  is  seen  that  the  effect  of  the  gate  at  each  point  in  channel  can  be 
modeled  by  introducing  capacitors  to  connect  the  channel  and  the  gate.  Therefore,  the  one¬ 
dimensional  equivalent  circuit  of  Fig.  5.2  can  be  extended  to  model  an  RG-FET.  If  the 
resistive  gate  is  assumed  to  be  ideal,  i.e.,  it  provides  the  linear  potential  drop  which  is 
independent  of  the  charge  oscillation  in  the  channel,  the  gate  can  be  modeled  as  a  series  of 
ideal  voltage  sources  as  shown  in  Fig.  5.6(a).  When  this  assumption  is  removed,  the  gate 

I 

can  be  modeled  as  a  series  of  resistors,  and  the  equivalent  circuit  Is  shown  in  Fig.  5.6(b). 


5.4.2  Qualitative  Description  of  the  RG-FET  Operation  | 

We  would  like  to  use  the  equivalent  circuit  in  Fig.  5.6  to  discuss  the  effects  of  the 
resistive  gate  on  the  operation  of  the  RG-FET  devices. 

First  of  all,  an  ideal  gate  is  assumed,  i.e.,  the  gate  provides  a  linear  potential 
distribution  independent  of  the  charge  oscillation  in  the  channel.  In  this  case,  the  gate  can 
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be  mcxleled  as  a  series  of  ideal  voltage  sources  (Fig.5.6(a)),  sinking  or  sourcing  as  much 
current  as  need  to  be.  When  tox  is  small,  the  gate  is  placed  very  close  to  the  channel,  the 
linear  gate  potential  distribution  has  great  control  over  the  channel  potential  and  the  gate  is 
able  to  screen  the  charge  variation  in  the  channel.  Contiguous  domain  oscillation  is 
expected.  As  the  gate  is  moved  farther  away,  Cox  becomes  much  smaller  and  the  gate  is 
basically  isolated  from  the  channel.  The  source  and  drain  gradually  take  over  and 
eventually,  single  domain  oscillation  is  expected. 

For  a  more  realistic  gate,  where  the  gate  is  modeled  by  real  resistors  (Fig.  5.6(b)),  the 
gate  voltage  is  no  longer  constant  The  gate  voltage  at  each  grid  point  is  determined  by  the 
channel  condition  and  has  to  be  solved  iteratively  in  the  simulation,  although  the  total 
potential  drop  on  the  gate  is  fixed  at  Vgq.  Current  loop  equations  have  to  be  set  up  for  that 
purpose.  As  seen  in  Fig.  5.6(b),  as  the  gate  resistance  gets  bigger,  the  displacement  current 
flowing  into  each  (Coj.  Cox)  branch  gets  smaller.  As  the  result  of  this,  more  displacement 
current  flows  along  the  channel,  which  means  source  and  drain  control  of  the  channel  is 
bigger.  Single  domain  oscillation  is  more  likely  to  occur.  On  the  other  hand,  as  the  gate 
resistance  is  reduced,  more  displacement  current  would  flow  into  the  gate,  giving  the  gate 
more  control  over  the  channel.  As  long  as  the  displacement  current  is  smaller  than  the  DC 
current  supplied  by  Vgg.  the  gate  potential  can  be  maintained  approximately  linear.  In  this 
case,  the  contiguous  domain  oscillation  is  expected. 


5.4.3  Simulation  Procedure 

In  the  simulation,  (<t)j,  nij)  is  chosen  as  a  set  of  pertinent  variables  for  the  discretized 
element  j,  where  <t>j  is  the  electrostatic  potential  depicted  in  Fig.  5.5(d)  and  nij  is  the  electron 
density  per  unit  width  in  the  channel. 

For  thi  equivalent  circuit  model  shown  in  Fig.  5.6,  the  total  charge  on  all  capacitor 
plates  is  zero  for  each  node.  At  node  j,  the  charge  balance  is: 

*  Qj.i+ Qj  +  Qnj  +  Qoj  -  ®  (5.35) 

Assuming  the  substrate  is  semi-insulating,  the  depletion  width  due  to  the  charge  on  the 
constant  capacitor  Ck  is: 


Q.  =  W  Ax.C„  ( <>.  -  <|).  ,) 

J  Vi"  (5.36) 

W  is  the  width  of  the  channel,  and  Axj  =  Tepi  -  x<ij  is  the  width  of  the  undepleted  channel 
region.  The  charge  on  Cn  is  derived  in  a  similar  way  as  described  in  equation  (5. 17): 
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(5.38) 


QNj  =  Wq(n,.-N^T^Ay) 

The  total  charge  on  the  channel  depletion  capadter  Cd  is  the  total  depletion  charge: 

Qn.  =  WAyqN„x^. 

(5  33) 

Equation  (5.34)  is  used  to  calculate  x^: 

*d=  /— AVj 

V  (5.39) 

(Jcj  in  (5.38)  is  a  non-linear  function  of  channel  potential  4>j.  To  linearize  the  Poisson's 
equation,  the  charge  in  the  depletion  region  (Joj  at  any  instant  of  time  is  expressed  in  terms 
of  the  charge  at  the  previous  time  via  the  Taylor  series.  Therefore, 
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<|)jOl(l  is  the  channel  potential  at  the  previous  time  and  jqjjoW  is  the  depletion  width  under  the 
gate  evaluated  at  (Pjo***. 
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V2«I»j-Vg/V;) 


Substituting  equations  (5.36-38)  into  (5.35),  the  set  of  linearized  Poisson's  equations 
for  a  doped  channel  RG-FET  with  m  elements  is: 

^1  ~  *^init 


“j-ifi-i-  Vl^j>  ♦)  *  “jVi*  ‘•v;]"*  'l(n,i-NoT^^iy)  - 


m  Yi„t 


(5.42) 
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Where 

k-T 

“  ^Gi'*'  ^gg‘‘‘  ^  ^“n"^ 


(5.43a) 

(5.43b) 


are  the  boundary  conditions  for  the  Poisson's  equation  assuming  the  source  and  drain 
contacts  are  ohmic,  a  and  p  are  defined  as: 


a  =  CjjAXj 

Equation  (5.41)  can  be  written  in  vector  notation: 


(5.44a) 

(5.44b) 


Ax  =  B 

where  A  is  a  tridiagonal  matrix.  Define  7J  =  aj-i  +aj  +  pj, 


(5.45) 


(5.46a) 


(5.46b) 
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(5.46)  can  be  easily  solved  for  a  given  electron  density  distribution,  assuming  an  initial  ^ 
distribution.  At  any  instant  of  time,  the  ^  distribution  has  to  be  iterated  until  a  consistent  ({> 
distribution  is  obtained. 

Once  the  potendal  distribudon  is  found,  the  electron  flux  can  be  calculted  and  the  new 
electron  density  an  instant  later  can  be  calculated  by  solving  the  current  continuity  equadon. 

The  boundary  condidons  for  the  electron  densides  per  unit  width  are: 


where 


nil  *  ^  ’^epi*  ’^dl^ 

n^  =  AyNn(T^..x^) 


(5.47a) 

(5.47b) 

(5.48a) 

(5.48b) 


AVdi  and  AVdm  are  obtained  by  subsdtudng  Vg=Vci  and  Vg=Voi+Vgg  into  (5.36), 
rcspecdvely. 

The  source  code  of  the  simuladon  program  is  listed  in  Appendix  E. 


5.4.4  Limitadons  of  the  RG-FET  Equivalent  Circuit  Model 

Examinadon  of  (5.35)  -  (5.38)  reveals  that  if  <!)j.i  »  ())j  (i.e.,  ([)  is  a  slowly-varying 
potendal  in  the  channel,  such  that  the  gradual  channel  approximadon  is  appropriate),  then 
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iL.  N_,Ay  ( T  X ..  ) 

“j  D  qji  dj'  (5.49) 

This  means  that  the  electrons  fill  the  undepleted  portion  of  the  channel  to  equal  the 
donor  density  Nd.  If  ^  <|>j  (i.e.,  changes  significantly  between  adjacent  grid  points), 
then  the  electron  density  in  the  undepleted  portion  of  the  channel  is 


n,j  =  ( V  (AXj.,+Ax.)  (t>.  +  Ax.<t)j^j) 


(5.50) 


n  - 

"  (T_,-  X  ..)Ay 

‘‘J  ^  (5.51) 

does  not  equal  to  Nd  and  the  channel  is  not  neutraL  Then  the  channel  cannot  be  reprssented 
by  a  single  potential,  as  in  Fig.  j.5,  and  the  lumped  equivalent  circuit  cannot  accurately 
represent  conditions  inside  the  device. 

Therefore,  we  have  to  be  aware  of  the  assumption  made  in  calculating  the  charge  Qd  on 
the  depletion  capacitor  Cd-  If  the  gradual  channel  approximation  is  valid,  i.e.,  the  gate  is 
dominant,  the  calculation  of  Qoj  is  accurate.  On  the  other  hand,  if  the  source  and  drain  are 
dominant,  the  contribution  due  to  Qd  is  significandy  smaller  than  the  Qj  terms  in  (5.35), 
and  the  inaccurance  introduced  in  calculating  Qd  would  not  affect  the  calculation.  When  the 
Qd  term  and  the  Qj  terms  are  comparable  in  magnitude,  we  might  anticipate  some  error. 


5.4.5  Simulation  Results 

The  objective  of  the  equivalent  circuit  model  is  to  simulate  both  single  domain 
oscillation  and  the  contiguous  domain  oscillation.  In  the  single  domain  mode,  the  image 
charge  of  the  domain  resides  on  sou";e  and  drain,  while  in  the  contiguous  domain  mode, 
the  image  charge  of  the  domain  resides  on  the  gate.  In  other  words,  the  channel  potential  is 
controlled  by  the  source  and  drain  in  the  single  domain  mode,  and  controlled  by  the  gate  in 
the  contiguous  domain  mode. 

The  equivalent  circuit  model  for  an  RG-FET  with  ideal  gate  can  be  used  to  examine  the 
effect  of  the  oxide  thickness  on  the  oscillation  mode.  The  thicker  the  gate  oxide  is,  the 
further  away  the  gate  is  from  the  channel,  the  greater  control  of  the  channel  potential  the 
source  and  drain  have,  and  more  likely  the  single  domain  mode  oscillation  will  occur.  On 
the  other  hand,  as  the  gate  oxide  thickness  is  reduced,  the  contiguous  domain  oscillation 
will  occur. 


mwi  iu»ii'yn!,iiijiiiijiiij|i^^ 
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The  following  simulations  arc  done  for  an  RG-FET  device  of  IC^m  long  and  24m 
wide.  The  channel  depth  is  1300A  and  the  doping  density  is  1.2xlO*'^ctn-3.  The  difference 
between  the  effective  metal  workfunction  and  the  electron  affinity  of  the  semiconductor 
-  X’ •*  0.5eV  and  the  bandgap  energy  is  1.36c V. 

5.4.5. 1  Gate  Dominated  Region  (tox  «  0 ) 

When  the  gate  oxide  is  aero,  the  device  is  an  RG-MESTET.  For  a  lOum  long  device,  a 
gate-to-gate  voltage  Vgq  of  8V  has  to  be  applied  to  bias  the  channel  into  the  negative 
mobility  region.  Fig.  5.7  shows  the  surface  electron  density  n,*nv*Ax  distribution  along 
the  channel  at  various  instants  of  time.  A  uniform  electron  density  is  assumed  initially 
except  in  the  middle  of  the  channel  where  a  dipole  domain  is  inmxluccd.  Very  quickly, 
random  oscillation  of  electron  density  occurs  in  the  channel.  This  oscillation  travels  towani 
the  drain  and  as  time  evolves,  the  oscillation  tends  to  stablize  (Fig.  5, 7(b),  (c)  and  (d)). 
After  approximately  one  device  transit  time  t  -  Uvavg  •  83ps  (vavg  is  typically  1.2x10^ 
cm/s),  the  oscillation  isstablizes,  as  seen  from  the  flux  vs.  time  plot  taken  in  the  middle  of 
the  channel  (Fig.  5.8).  The  quasi-steady  state  electron  density,  electric  field  and  potenti.il 
distribution  along  the  channel  at  time»100ps  is  shown  in  Fig.  5.9.  In  Fig.  5.9(.i), 
contiguous  domains  of  surface  electron  density  art  observed.  The  traveling  vekxiiy  of  the 
domain  can  be  found  by  tracing  the  sanx  domain  at  two  different  times,  which  is  1.3x10^ 
cm/s  for  this  case.  The  oscillation  frequency  is  lOOGHz,  as  determined  from  Fig.  5.8,  As 
the  electron  density  oscillates,  the  electric  field  and  the  potential  also  oscillate,  as  shown  in 
Fig.  5.9(b)  and  (c).  . 

The  oscillation  frequency  and  peak-to-peak  amplitude  depend  on  the  biases  when 
tox"0A.  Vso  "  Vdo  values  varying  from  0  to  I.2V  have  been  investig.itcd,  while  Voo  is 
held  at  8V.  In  Fi^.  5.10,  the  simulation  results  from  the  equivalent  circuit  mcKlcI  (black 

i 

dots)  are  superposed  to  the  simulation  results  for  an  RG-MESFET  from  the  TOF  program 
(Fig.  3.9),  and  the  two  results  are  very  similar,  except  at  higher  l/n,X(i  values. 
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Figure  5.7  Tlie  time  evolution  of  electron  density  in  the  channel  of  an  RG-MESFET. 
biased  at  Vgg*8V,  Vgi®0V,  Vsg*0-8V,  Vdo^O.SV. 
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Simulation  results  ((a)  ns,  (b)  ^y,  and  (c)  (j>  along  the  channel)  at  t^lOOps 
for  an  RG-MESFET  (tox=0)  using  the  equivalent  circuit  model.  Bias 
condition;  Vgg=8V,  Vni^V,  Vsg^-8V,  Vix}=0.8V. 
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Figure  5.10  Comparison  of  simulation  results  using  the  equivalent  circuit  model  (big 
black  dots)  with  that  of  time-of-flight  simulation  (Fig.  3.9).  (a)  frequency, 
(b)  amplitude. 


S.4.S.2  Souice/Drain  Dominated  Region  (tox  »  L) 

On  the  other  hand,  when  tox  “  6OOOOOA,  the  gate  is  so  far  away  from  the  channel  that  it 
can  not  influence  the  channel  potential.  The  channel  potential  is  soley  determined  by  the 
source  and  the  drain.  Only  the  single  domain  oscillation  is  observed  in  this  case.  As  seen  in 
Fig.  5.1 1,  the  curves  for  tox  =  600(X)0A  are  very  similar  to  that  of  Gunn  diode  shown  in 
Fig.  5.3.  The  domain  travels  with  a  velocity  of  about  7-5x10®  cm/s. 
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Figure  5.11 
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Simulation  resulats  ((a)  n*,  (b)  §y,  and  (c)  6  along  the  channel)  at  different 
times  for  an  RG-MOSFET  (tox=600000A)  using  the  equivalent  circuit 
model.  Bias  condition:  V(3g=8V,  Vgi=OV,  Vsg=0-8V,  Vix}=0.8V. 
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5.4.5.3  Intermediate  tox  Region 

The  simulation  program,  when  nm  for  intermediate  tox  values  between  0  and  6OOOOOA, 
exibits  convergence  problems.  This  may  be  due  to  the  failure  of  the  gradual  channel 
approximation,  as  discussed  in  section  5.4.4.  If  the  gradual  channel  approximation  fails, 
the  simple  picture  of  the  channel  electrostatics  given  in  Fig.  S.S  breaks  down,  and  a  more 
sophisticated  modeling  approach  is  needed.  It  is  also  possible  that  the  convergence 
problems  may  simply  be  due  to  the  numerical  algorithms  employed  in  the  simulation. 
Further  investigation  of  the  intermediate  tox  region  is  beyond  scope  of  the  present  thesis, 
but  is  an  interesting  area  for  future  development 


5.5  Summary 

The  equivalent  circuit  model  successfully  simulates  the  RG-FET  devices  in  the  two 
extremes:  tox=0  and  tox*>^<  When  the  gate  is  placed  very  close  to  the  channel,  the  channel 
potential  is  mainly  determined  by  the  gate  volt<*<!e.  Contiguous  domain  oscillation  is 
obtained.  On  the  other  hand,  when  the  gate  is  placed  far  away  from  the  channel,  the 
channel  potential  is  mainly  determined  by  the  source  and  drain.  Single  domain  oscillation  is 
obtained.  This  general  behavior  is  consistent  \dth  the  t>  variable  experimental  results  and 
with  expectations  based  on  physical  intuition. 
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CHAPTER  6 

CONCLUSION  AND  RECOMMENDATION 


The  resistive  gate  FET  devices  are  unique  in  many  aspects.  By  applying  gate-to-gate 
voltage  on  the  resistive  gate,  high  and  unifam  electric  field  can  be  obtained.  With  the  right 
gate  sheet  resistance,  the  gate  can  screen  the  space  charge  in  the  channel  and  contiguous 
domain  oscillation  can  occur.  This  research  project  provides  the  unique  opportunity  to 
further  understand  the  resistive  gate  FET  devices  both  experimentally  and  theoretically. 

Experimentally,  the  fabricated  devices  have  been  DC  and  microwave  tested.  Microwave 
oscillations  are  observed,  but  only  in  the  single  domain  mode.  Preliminary  analysis 
indicates  that  the  gate  sheet  resistaiice  has  to  be  reduced  to  about  SkD/square  to  obtain 
contiguous  domain  oscillation.  Cermet  deposition  was  investigated  to  obtain  a  desirable 
resistive  film  for  future  fabrication  of  CDO  devices. 

Theoretically,  the  steady  state  simulations  of  the  long  channel  RG-FET  devices  with 
constant  mobility  have  been  performed.  The  results  indicate  that  the  operation  of  the  RG- 
FET  experiences  three  regions  for  positive  current  They  are  divided  by  the  uniform 
channel  condition  and  the  pinch-off  condition.  The  pre-uniform  region  is  solely  caused  by 
the  non-zero  Vgg-  The  post-uniform  region  and  the  pinch-off  regions  are  very  similar  to 
the  conventional  FETs  ohmic  and  pinch-off  regions. 

Transient  simulations  of  RG-FET  devices  have  also  been  performed.  By  combining  the 
transferred  electron  effect  in  GaAs  with  the  strong  gate  control  of  the  channel  potential,  a 
series  of  contiguous  electron  packets  can  be  generated,  as  indicated  by  the  one  dimensional 
computer  simulations.  High  frequency  oscillation  in  the  GHz  range  can  be  obtained.  The 
fiequency  can  be  tuned  by  changing  the  source  to  gatel  voltage. 

To  understand  the  single  domain  oscillation  observed  in  experimental  RG-MESFET 
devices,  an  equivalent  circuit  model  was  developed.  Two  modes  of  oscillation  are  obtained. 
When  the  gate  is  placed  very  close  to  the  channel,  the  channel  potential  is  mainly 
determined  by  the  gate  voltage,  and  contiguous  domain  oscillation  is  observed.  On  the 
other  hand,  when  the  gate  is  placed  far  away  from  the  channel,  the  channel  potential  is 
mainly  determined  by  the  source  and  drain,  and  single  domain  oscillation  is  observed. 
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There  are  still  many  questions  to  be  answered.  The  future  work  should  consist  of  two 
parts.  First  of  all,  new  devices  should  be  fabricated  with  reduced  gate  sheet  resistance.  The 
new  devices  should  eliminate  the  resistive  gate  overlapping  the  source  and  drain  region,  so 
to  eliminate  the  coupling  between  Vgg  and  Vsg-  More  test  structures  should  be  put  on  the 
die  to  better  monitor  the  fabrication  process. 

The  theoretical  work  is  essential  to  further  understanding  of  the  device.  Better 
approximation  of  the  charge  on  the  depletion  capacitor  Cd  should  be  explored  so  the  model 
can  properly  represent  the  channel  when  ny  ^  Nid  and  simulate  the  intermediate  tox  values. 
Once  this  step  is  tinished,  the  equivalent  circuit  model  should  be  extended  to  include  the 
resistive  gate  effect  The  finite  resistance  of  the  resistive  gate  tends  to  alter  the  gate  field, 
thus  alter  the  electric  field  in  the  channel.  Upon  Hnishing  the  resistive  gate  part,  the 
theoretical  work  would  come  to  a  logical  conclusion. 

A  full  2D  simulation  program  might  be  needed  to  correctly  model  the  RG-FET  device. 
However,  the  commercially  available  2D  device  simulator  PISCES  also  has  convergence 
problems  when  dealing  with  the  nonlinear  velocity-field  curve  of  GaAs. 
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Appendix  A 

General  Channel  Potential  Expressions 
for  MOS  Structure  Using  Green's  Function 


This  appendix  calculates  the  channel  potential  of  MOSFET  using  Green's  function[20]. 
Upon  appropriate  approximation,  it  gives  the  potendal  in  the  channel  as  well  as  the  non¬ 
local  potential  correction. 


Figure  A-1  MOS  system.  The  potential  at  any  point  in  the  semiconductor  (x,y)  due  to  a 
line  charge  parallel  to  the  z  axis  at  (x',y')  is  calculated  using  the  method  of 
image. 


Consider  a  MOS  structure  as  shown  in  Fig.  A-1.  The  potential  at  any  point  (x,y)  in  the 
semiconductor  region  (with  dielectric  constant  £$)  due  to  a  linear  charge  of  unit  strength, 
i.e.,'  IC/cm,  at  a  point  (x',y')  parallel  to  the  z-axis  can  be  calculated  by  the  method  of 
images: 


G(x,y,x',y’) 


-1 


ln((y-y*)^+  (x+x')^  +  ln[(y-y-)^+  (x-x’)^  ■ 


"OX 


in=l 


1 

^  6.-6  ' 
^  ^ox 

2a  1 

Vn-1 


ln[(y-y’)^+(x+x’+2mtox)^] 


(A-1) 


where  tox  is  the  oxide  thickness.  This  potential  function  is  the  Green  function  solution  of 
the  two-dimensional  Poisson  equation  of  the  structure.  The  metal  plane  is  assumed  at 


ground  potential.  Since  the  boundary  conditions  at  both  the  insulator-semiconductor 
interface  and  insulator-metal  interface  should  be  satisfied,  an  infinite  series  of  image 
charges  are  'squired  to  calculate  the  potential  function  in  each  regio". 

Define: 


c  ~ 

0  4jie,(e,+eox) 

tl-l 

®ox 

2 

(A-2a) 

7c(e,+eox) 

(A-2b) 

am  =  x-t-x’+2mt„. 

(A-2c) 

Eqn  (A-1)  can  be  simplified  as: 

G(x,y,x*,y*)  =  ln[(y-y')*+  (x-x')^  +  Cn,ln[(y-y)^+  a^] 

*  (A-3) 

When  y'  approaches  infinity,  the  potendal  should  be  zero.  In  (|A-3),  all  nature  log  functions 

approach  infinity  at  the  same  rate.  The  sum  of  the  coefficients  can  be  proven  to  be  zero.  i.e: 

— 

y  c„=-^ 


m>0 


47:e, 


(A-4) 


Physically,  it  means  that  the  sum  of  the  real  charge  and  all  the  image  charges  is  zero. 


Appendix  B 

Electrostatics  of  the  RG-MOSFET 
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Assuming  a  MOSFET  structure  ./ith  oxide  thickness  tox.  substrate  (p-type)  doping  Na. 
The  surface  potential  at  any  surface  point  (0,y)  in  a  MOSFET  is  simply  the  sum  of  ideal 
gate  voltage  Vq',  the  potential  due  to  the  inversion  charge  and  the  potential  due  to  the  bulk 
charge  (refer  to  Fig.  A- 1): 

• 

Vs(y)  =  V^Cy)  +jG(0,y,0,y')(-qn,(y‘))dy’  +Jdy‘jG(0,y,x’,y’)(-qN^)dx- 

-  -  0  (B-1) 

where  X(i(y')  =  depletion  layer  thickness  defined  by  the  implicit  relationship: 


Vs(Xd.y)  =  V^Cy)  -^jG(x^,y,0,y’)(-qn,(y’))  dy 


+  Jdy’jG(x^,y,x’,y>(-qN^)  dx*  =  0 


—  0 


If  nsCy*)  and  xd(y')  are  taken  to  be  constants,  ns(y’)  in  second  term  can  be  taken  out  of 
the  integral  and  the  integration  with  respect  to  x'  and  y'  can  be  interchanged  in  thin!  term.  It 
can  be  shown  that  (Appendix  D) 

M 

f  G(x,y,x’,y*)  dy*  =  |-  +  x  >  x’ 

®  ^  X  <  x’ 

*  (B-3) 

Eqn  (B-1)  and  (B-2)  become 

V,  =  VQ(y).3^..3^ 

'-ox  '-ox 


Solving  xd  from  Eqn  (P-5),  Eqn  (B-4)  becomes: 

¥.0-;  -  Vi(y)  -  5^  +  V„-  /  vJ+2V„(V;.(y) 


qns(y) 


which  is  exactly  the  same  as  the  one-dimensional  solution  of  the  Poisson  equation  using  the 
depletion  approximation  for  the  MOS  smicture(9].  Here 


Since  the  surface  charge  in  the  channel  of  a  real  MOSFET  is  not  a  constant,  the  non¬ 
local  potential  correction  is: 

M 

AV,(y)  =  -qjG(0,y,0.y‘)  ( n,(y’)  -  n,(y) )  dy’ 

=  q  ff  in(yy‘^ 

2n(e,+eox)JI,  e.+Sox  j 

( n,(y*)  -  n,(y) )  dy’ 

Hiis  is  accurate  so  long  as  the  variation  in  the  depletion  width  is  neglected  (the 
depletion  width  xd  is  pinned  once  the  MOSFET  enters  the  inversion  region). 


Appendix  C 

Electrostatics  of  the  RG-MESFET 
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The  MES  structure  has  no  oxide  underneath  the  gate  and  there  is  a  doped  channel  with 
doping  density  of  Nd  and  channel  depth  Tepi-  Before  the  channel  is  pinched  off,  the  cross 
section  of  this  RG-MESFET  can  be  viewed  as  a  Schottky  diode  in  series  with  an  n'^p 
junction  diode  (as  shown  in  Fig.  C-1).  The  electrons  fill  up  in  between  the  two  diodes  with 
a  volume  density  equal  to  the  channel  doping  Nq. 


Figure  C-1  MES  system.  The  potential  at  any  point  in  the  channel  (x,y)  due  to  a  line 
charge  parallel  to  the  z  axis  at  (x'.y')  is  calculated  using  the  method  of 
image. 


The  Green  function  for  a  MES  structure  can  be  obtained  by  setting  es=eox  and  t  ox=0  in 
Eqn  (A-1)  of  Appendix  A: 

G(x,y,x’,y‘)  =-r^[  InCCy-y’j^+Cx-x*)^  -  ln[(y-y’)^+(x+x')^) 

(C-l) 

The  sum  term  in  Eqn  (A-1)  is  eliminated  for  there  is  only  on  image  for  each  line  charge. 

The  channel  potential  is  the  sum  of  the  ideal  gate  voltage  V’o  and  potential  due  to  the 
depleted  channel  and  substrate; 


TT 
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ijititaa<wa»^*«(ife>i A  liwin'iiiiBi^ 
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«  *<(y  ) 

Vti,(y)  -  VQ(y)  ♦  jdy  jccx'.y’.x.y)  qN^d 


—  0 


T  ♦«-(>•) 


■  Jdy’JoCx’.y’.x^)  qN^dx*  ♦  Jdy’jG(x\y\x,y)(-qN^)dx‘ 


—  T 


where 


(C-2) 


(C-3) 


<^Hn  i*  the  barrier  heighc  of  the  Schottky  duxle.  and  xd(y')  the  depiction  layer  thickness  of 
the  Schottky  diode.  Xd2(y')  and  Xp(y*)  are  the  depletion  width  of  the  channcl-iubst  atc 
junction  diode. 

If  Xd(y'),  x<i2(y')  and  Xpfy*)  are  constants  independent  of  y’.  i.e.  charge  in  the  channel 
are  uniform  across  the  channel,  the  Integration  with  respect  to  y'  and  x’  can  be 
interchanged.  Using  lujn  (B-3),  Hqn  (C-2))  becomes; 

v.>,(y)  -  v^fy)  +  ,  J  ♦ 

-V„(y)^^B.x] 

(C-4) 

The  last  two  *^nns  cancel  each  other  due  to  the  fact  that  tlie  negative  charge  in  die  substrate 
must  be  equal  to  the  positive  charge  in  the  channel  for  the  channel-substrate  junction- 

W  *<t2>"  Vp 


C,  P 


Surface  electron  density  is  ilrfined  to  be: 

n.(y)  -  Xj(y).x,^2^y)) 

By  solving  Fqn  (C-3),  (C  6)  together  with 

«  S<y  > 


(C  5) 
(C6) 


■M  C.(x'.y'.T,^*.Xp,y)  qN„dx' 


+  Jdy  J  G(x'.y'.T^p.4-Xp.y)  (  qN^)  dx’  »  0 


--  T 


xj  is  found: 


«A— *, 


Nn+N.  2  2  N„+N.  2£.N.  ^ 

(■  p  -  A/t  .  P  ^  t’  4-  *  A  v* 

'  ISI  ^  ‘  M  ‘  2  0 
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where 


t‘-T 


•P*  Nr 


For  a  semi-insulating  substrate,  Na=0,  the  channel  potential  is: 

Vjy)“VG(y)  +  ^(T  --1^) 

I  D  (C-10) 

The  non-local  potential  correction  is  caused  by  non-uniform  charge  in  the  channel  (position 
dependent  x<i(y).  Xd2(y)  and  Xp(y)).  Still  assuming  that  the  doping  in  the  substrate  is  zero, 
i  e.,  x<i2(y)*=Tept,  the  non-local  potential  will  only  arise  from  the  variation  in  x<i(y). 
MathcmaticJiIly, 

- 

Av^^(x.y)  «  qNpJdy'  J  G(x.y,x’,y’)  dx’ 

'  (C-ll) 

The  integration  with  respect  to  x’  and  y'  can  no  loner  be  interchanged.  The  integration  with 
res;>cct  to  x’  is: 


J  G(x,y.x’,y  )  dx’  -  [  (x-x^(y))  In  ((y-y’)V  (x-Xj(y))^ 

*/y)  * 

-  (X-Xd^y’))  In  ((y-y’)^+  (X-X^(y’))^ 

-  (x+x^(y’))  In  ((y-y’)^+  (x+Xj(y’))^ 
(x+Xj(y))  In  ((y-y’)^+  (x+Xj(y))^ 

+  2  ly-y  I  (  un  -  tan  * 
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Proof  of  Equation  (B-3) 
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Equation  (B-3)  in  Appendix  B  is: 

m 

jG(x.y.x*y)dy-«^  +  ^ 


x>x’ 


where  G(x,y,x’,y’)  is  expressed  in  Eqn  (A-3)  of  Appendix  A,  repeated  here: 

G(x.y.x’.y’)  =  ^  ln[(y-/)^  (x-x’)^  +  £  Q„ln[(y-y*)\  a^ ] 

*  m-O 

with 
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Cm «  j  6|+£ox J  “  2  ^ 
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am  -  x  +  X’  +2mtox 

First  let’s  calculate  a  simple  integratin: 

«  y 


J  ln[(y-y*)^+  b^]  dy  »  Jln[(y-y*)^+b^]dy’  +J  ln[(y.y’)^-t-b^]dy' 

m 

» ij  ln[u^+b^  du 
^  2u  ln(uVb^  -4u  +4lbl  tan'* 


Apply  Eqn  (D-4)  to  (D-2): 


(D-3a) 


(D-3b) 

(D-3c) 


jG(x,y.x*.y')dy’  -  2u  Infu^+fx-x’)^  -4u  +4  Ix-x'i  tan-»^|^ 

+  ]£  [  2u  in(u2+a2  )  .4u  +4am  tan-*^J~ 

m-o  ^  (J5.5) 

Note  that  when  y'  approaches  infinity,  ail  mature  log  functions  approach  infinity  at  the 
same  rate.  Using  Eqn  (A-4)  of  Appendix  A,  (D-5)  becomes: 


m  ^ 
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*  m»0 


In  the  case  of  x  >  x' 


Since 
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Therefore, 


f  G(x,y,x’.y’)dy’  +  ^  =  +  ^ 

J  V-OC 


(P-6) 

(P-7) 

(D-8) 


(D-9) 


Case  for  X  <  x'  can  be  obtained  in  a  similar  way. 


